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Summary 


An analytical, parametric study of the attenuation of bending boundary layers or edge effects in 
balanced and unbalanced, symmetrically and unsymmetrically laminated thin cylindrical shells is pre- 
sented for nine contemporary material systems. The analysis is based on the linear Sanders-Koiter shell 
equations and specializations to the Love-Kirchhoff shell equations and Donnell’s equations are 
included. Two nondimensional parameters are identified that characterize and quantify the effects of 
laminate orthotropy and laminate anisotropy on the bending boundary-layer decay length in a very 
general and encompassing manner. 

A substantial number of structural design technology results are presented for a wide range of 
laminated-composite cylinders. For all laminate constructions considered, the results show that the 
differences between results that were obtained with the Sanders-Koiter shell equations, the Love- 
Kirchhoff shell equations, and Donnell’s equations are negligible. The results also show that the effect 
of anisotropy in the fomi of coupling between pure bending and twisting has a negligible effect on the 
size of the bending boundary-layer decay length of the balanced, symmetrically laminated cylinders 
considered. Moreover, the results show that coupling between the various types of shell anisotropies has 
a negligible effect on the calculation of the bending boundary-layer decay length in most cases. The 
results also show that, in some cases, neglecting the shell anisotropy results in underestimating the 
bending boundary-layer decay length and, in other cases, results in an overestimation. An example 
problem is included in an appendix that demonstrates how to perfomi the calculations that were used to 
generate the results of the present study. 

Symbols 

a 12> a 22’ a 26’ B 2l inverted stiffness expressions defined in appendix C 

a 12’ a 22’ a 26’ B 21 modified inverted stiffness expressions defined in appendix B 

nondimensional anisotropy parameters 

An, A i2’ A i6 , laminate membrane stiffnesses 
A 22’ A 26’ A 66 

a 16’ A 26’ A 66 modified laminate stiffnesses defined in appendix B 

Bn, Bi 2 , B i6 , laminate membrane-bending coupling stiffnesses 
B 22’ B 26’ B 66 

b 16’ B 26 modified laminate stiffnesses defined in appendix B 

first-order correction factor for anisotropy parameter 

d, d° attenuation or decay lengths for which anisotropy is included and neglected, 

respectively 

Du, Di 2 > D 16 , laminate bending stiffnesses 
D 22> D 26> D 66 

e 0 


stiffness coefficients defined in appendixes C and B, respectively 
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major and minor principal lamina moduli, respectively 
lamina shear modulus 

cylinder wall thickness and length, respectively 

axial and circumferential bending and twisting stress resultants, respectively 
axial, circumferential, and shear membrane stress resultants, respectively 
nondimensional orthotropy parameter 

loading function appearing in bending boundary-layer differential equation 
constant coefficients of bending boundary-layer differential equation 
axial and circumferential transverse- shear stress resultants, respectively 
cylinder radius 

modified shear stress resultant defined in appendix B 
strain-energy density 
radial-displacement component 
axial coordinate of cylinder 
attenuation-length tolerance parameter 

axial, circumferential, and shear membrane strains, respectively 
circumferential, angular coordinate 

axial and circumferential bending and twisting strains, respectively 
constant defining different shell theories 
lamina major Poisson’s ratio 
lamina fiber angle (see fig. 4) 


Introduction 


The term, “bending boundary layer,” refers to localized zones of bending stresses and deformations 
that appear in practically every type of thin-shell structure (ref. 1). Bending boundary layers are caused 
by edge support conditions; by localized mechanical loads, heating, or cooling; and by abrupt changes 
in stiffness, such as those caused by a cutout, a crack, or a stiffener. All these effects may be real con- 
cerns in a given preliminary design for an advanced aerospace vehicle made of laminated-composite 
materials. Thus, it is useful to have nondimensional parameters that characterize the effects of shell 
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geometry and laminate orthotropy and anisotropy on the extent of bending boundary layers that can be 
used to help guide the development of a design. For example, an optimal design for a pressure vessel 
might be one that exploits the membrane load-carrying action of a shell and minimizes zones of local 
bending stresses. Meaningful nondimensional parameters could be used at the preliminary design stage 
to identify families of laminates and material systems that exhibit relatively small bending boundary 
layers. Moreover, a meaningful estimate of the size of a bending boundary layer in a shell is very useful 
for determining an adequate first-approximation finite-element model for a complex shell structure. 
Without a proper understanding of the extent of a bending boundary layer, it is possible to have a finite- 
element model that could miss a significant part of the structural response in a region where failures are 
often initiated by high interlaminar stresses. Furthermore, a priori knowledge of the extent of bending 
boundary layers is useful in determining the instrumentation locations in structural verification tests and 
in material characterization tests (ref. 2). In addition, knowledge of how laminate construction affects 
the extent of a bending boundary layer is useful for understanding how nonlinear prebuckling deforma- 
tions affect the buckling behavior of cylindrical shells. 


Studies of the behavior of axisymmetric, bending boundary layers in right-circular, cylindrical shell 
structures made of orthotropic or anisotropic materials and with finite length have been presented, to at 
least some extent, in references 1 through 13. In the discussion that follows, reference is made to unbal- 
anced and balanced laminates that are either symmetrically or unsymmetrically laminated. Herein, the 
temi unbalanced laminate is used to indicate that coupling between pure extension or contraction and 
shearing is present in a laminate. The term unsymmetric laminate is used to indicate coupling between 
any of the components of bending action with any of the components of membrane action. A fully 
anisotropic laminate would include both these types of anisotropy, in addition to the anisotropy that is 
manifested by coupling between pure bending and twisting action that is sometimes exhibited by 
balanced, symmetric laminates. 


In reference 1, an analysis is presented, and an expression for the attenuation or decay length of the 
bending boundary layer for a specially orthotropic cylinder that is subjected to edge loads, internal 
pressure, and heating is given. These equations, and the accompanying results, are based on the linear 
Love-Kirchhoff shell equations. In reference 3, an analytical solution that is based on DonnelFs simpli- 
fications to the linear Love-Kirchhoff shell equations is given for fully anisotropic cylinders that are 
subjected to lateral pressure and edge loads. Results that show the effect of laminate anisotropy on the 
edge moment are presented for a clamped two-ply shell that is subjected to internal pressure. In 
addition, a discussion is presented that suggests that solutions based on DonnelFs equations should be 
accurate for laminates that are not highly anisotropic. An analytical study of bending boundary layers in 
unbalanced, symmetrically laminated cylinders, which is also based on DonnelFs equations, is 
presented in reference 2. The aim of this study was to determine a suitable gage section in a laminated- 
composite tube that is to be used for a material characterization test. Results are presented for 
unidirectional, helical-wound tubes. 


An analytical solution for bending boundary layers in unbalanced, symmetrically laminated and 
balanced, unsymmetrically laminated circular cylindrical shells that are subjected to internal pressure 
and themial loads is presented in reference 4. The solution is also based on DonnelFs linear equations, 
and numerical results are presented for filament-wound cylinders made of heat-treated carbon-carbon 
material. A study that focuses mainly on prebuckling deformations, with bending boundary layers, in 
homogeneous, orthotropic and unsymmetrically laminated cross-ply cylinders that are subjected to 
axial-compression loads and lateral pressure loads is presented in reference 5. The effects of the bend- 
ing boundary layers on the buckling response are examined for several laminate constructions, but the 
general effects of the laminate construction on the extent of the boundary layers are not discussed. 
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A pair of complex conjugate, fourth-order equations that are based on Flugge’s corresponding 
equations (ref. 14), which can be solved in closed form, are derived for specially orthotropic, circular 
cylindrical shells in reference 6. Moreover, eigenfunction solutions are presented that include the solu- 
tion for the axisymmetric bending boundary layer; several simplified equations are presented and their 
relative accuracy is analyzed. In reference 7, a study of bending boundary layers in transversely isotro- 
pic circular cylindrical shells is presented. This study examines the attenuation characteristics of bend- 
ing boundary layers by applying an asymptotic method to the linear, three-dimensional elasticity 
equations and presents order-of-magnitude estimates for the stresses and displacements for a wide range 
of ratios of the two principal elastic moduli. In reference 8, an analytical solution for an unbalanced, 
unsymmetrically laminated circular cylindrical shell that is subjected to internal pressure is presented 
that is based on a variant of the Love-Kirchhoff shell theory, which uses an expression for the change in 
surface twist that was given by Timoshenko. Numerical results that demonstrate the coupling effects of 
the shell anisotropies are also presented for a two-ply shell. 

The bending boundary layers of an unbalanced, unsymmetrically laminated circular cylindrical 
shell that is subjected to axial compression, torsion, or thermal loading are investigated in reference 9. 
Results are also presented that demonstrate the coupling effects of the shell anisotropies. In addition, 
results are presented for two more conventional unsymmetric laminates and a typical quasi-isotropic 
laminate. In references 10 and 11, bending boundary layers are also examined for balanced, 
symmetrically laminated and balanced, unsymmetrically laminated cylindrical shells, in the context of 
nonlinear prebuckling deformations that occur as a result of compression and thermal loads. In 
particular, the effects of laminate stacking sequence on the extent and character of the bending 
boundary layers are presented for two groups of three similar laminates. Two of the laminates are 
unsymmetric. In reference 12, a linear analysis is presented that focuses mainly on balanced, symmetri- 
cally laminated cylinders, and an expression is given for the length of the bending boundary layers near 
the cylinder ends that is based on the Love-Kirchhoff shell equations. 

Most recently, Goldenveizer’s static-geometric duality principle (ref. 15) has been used in refer- 
ence 13 to reduce the Sanders-Koiter equations (refs. 16 and 17) for fully anisotropic, right-circular 
cylindrical shells to two coupled, fourth-order equations that use a stress and a curvature function as the 
unknown, primary field variables. The reduction is done by adding certain negligibly small temis to the 
stress-strain relations, which are intrinsically in error because they must be established experimentally. 
The approach demonstrates how the static-geometric duality principle can be used to reduce greatly the 
amount of algebra needed to obtain results. Eigenfunction solutions are also presented for specially 
ortho tropic cylinders that are in agreement with corresponding results presented in reference 6. More- 
over, asymptotic formulas that can be used to determine axisymmetric bending boundary-layer attenua- 
tion lengths and the decay of other unsymmetric, self-equilibrated edge loads are given. 

With the exception of reference 13, explicit expressions for estimating the size of axisymmetric 
bending boundary layers in fully anisotropic, right-circular cylinders are not found in the literature. 
Moreover, there appear to be even fewer results for laminated-composite shells made of contemporary 
material systems and essentially no substantial parametric studies. The present paper focuses on devel- 
oping meaningful estimates of attenuation lengths of bending boundary layers in balanced and 
unbalanced, symmetrically and unsymmetrically laminated circular cylinders. The analysis is based on 
the linear Sanders-Koiter shell equations and contains the Love-Kirchhoff shell equations (ref. 1) and 
Donnell’s equations (ref. 1) as special cases, and is somewhat similar to the analyses presented by 
Reuter (ref. 4) and Chaudhuri, Balaraman, and Kunukkasseril (ref. 8). With these equations, explicit 
expressions are obtained and nondimensional parameters are presented that characterize the effects of 
cylinder geometry and laminate construction on the size of a bending boundary layer in a very general 
manner. In particular, generic design curves are presented that use the nondimensional parameters to 
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show the effects of laminate orthotropy and anisotropy on the attenuation length in a concise and 
encompassing manner. In addition, values of these parameters are presented for a very wide range of 
ortho tropic and anisotropic laminate constructions. Also, differences in the results that were obtained in 
the present study by using the Sanders-Koiter shell equations, the Love-Kirchhoff shell equations, and 
Donnell’s equations are discussed. Key details of the analysis presented are elaborated upon in 
appendixes A-C. An example problem is included in appendix D that demonstrates how to perfomi the 
calculations that were used to generate the results of the present study. 


Analysis 

The ordinary differential equation that governs the axisymmetric bending behavior of a right- 
circular cylinder that is subjected to edge loads or displacements and surface tractions is obtained by 
first specializing the linear Sanders-Koiter shell equations that are given in appendix A for axial sym- 
metry. For the equations presented herein, x and 0 denote the axial and circumferential coordinates of 
a right-circular cylinder, respectively, and the specialization to axial symmetry is obtained by eliminat- 
ing all temis in the equations that are differentiated with respect to the circumferential coordinate 0. The 
resulting set of equations for axisymmetric behavior is given in appendix B. The ordinary differential 
equation that governs the axisymmetric bending behavior of a right-circular cylinder that is subjected to 
edge loads or displacements and surface tractions is derived in appendix B and is given by 

d 4 w d 2 w „ _ . 

—7 +4S — +4Qw =P(x) ( 1 ) 

dx 4 dx 2 

where S, Q, and P(x) are defined in appendix B by equations (B55), (B56), and (B57), respectively, 
and w(x) is the radial deflection that is positive valued when outward. The coefficients of equation (1) 
depend on the subscripted A, B, and D constitutive terms of classical Love-Kirchhoff-type laminated 
shell theory (e.g., see ref. 18, pp. 190-202) and the radius of the cylinder middle surface R. 

To detemiine the specific fomi of the solution to equation (1), it is useful to examine the positive- 
definiteness conditions on the strain-energy density function. The strain-energy density function for this 
problem is given by 


I'll - N x £° + N e 80 + N x0 ^0 + M x k° + M x 0 K °0 ( 2 ) 

where N x , N 0 , and N x0 are the membrane stress resultants, M x and M x0 are the bending stress 
resultants, 8 X , 8 0 , and t £ 0 are the middle-surface membrane strains, and K x and K x0 are the middle- 
surface bending strains. By using equations (B22), (B23), and (B28), the strain-energy density function 
is expressed as 


'I'll = N x e° + N e 80 + Ty ^ 0 + M x k° 


( 3 ) 
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The strain-energy density is expressed in terms of the strains and constitutive terms by using the consti- 
tutive equation given by equation (B29); that is, 
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The stiffness terms in equation (4) that have overbars are defined by equations (B31)-(B35) and are 
functions of the shell wall thickness-to-radius parameter, — . By enforcing positive definiteness of the 
strain-energy density function (e.g., see ref. 19), the requirement that the diagonal terms A 11? A 2 2 , 
A 66 , and D n be positive valued is obtained. Moreover, the following determinants are positive valued: 


A 11 A 12 
A 12 A 22 


- A 11 A 22~ A 12 > ® 


(5) 
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A 66 - A 11 A 26 “ A 22 A 16 + 2A 12 A 16 A 26 > 0 


( 6 ) 


Likewise, positive definiteness of the strain-energy density function also requires that the determinant 
of the constimtive matrix in equation (4) be positive valued. Moreover, by rearranging the strain-energy 
density function into the form 
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the following additional positive-definiteness condition is obtained: 


A 

A 


11 

16 


46 


l 66 


-A 11 A 66 - A 16 > ^ 


(7) 


( 8 ) 


6 



The homogeneous solution for equation (1) involves the square root of the quantity Q-S“ 
using equations (B55) and (B56), this quantity is given by 


By 


Q-S 2 = 


4C!C 3 -Ci 

16C 2 


(9) 


2 

Substituting equations (B41)-(B43) into equation (9) and simplifying, the quantity Q- S~ is found to 
be given by 
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( 10 ) 


2 

It follows logically that Q- S w >0 because the positive-definiteness of the strain-energy density func- 

2 

tion requires that the detemiinant in equation (10) be positive valued. Moreover, Q-S w >0 implies 

C 3 

that Q > 0, and Q > 0 implies that >0. Equations (6), (8), and (B43) indicate that C 3 >0. Thus, 
C 3 Cl 

— >0 yields the condition that C 1 -D L1 e>0 (see eqs. (B49) and (B50)). Because >0; e> 0. 

Ci 

Because Q is always positive, it is convenient to introduce the expression 


t 2 =q= — 

4R 2 a 22 D U e 
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such that T - S > 0 , and to express equation (1) as 


d 4 w _ d 2 w ™ x 

j + 4S — - +4Tw = P(x) 

dx 4 dx 2 


( 12 ) 


Equation (12) is a linear, fourth-order, nonhomogeneous ordinary differential equation with con- 
stant coefficients. The characteristic equation of equation (12) is given by 


X 4 + 4S ^ 2 + 4T 2 = 0 


(13) 


n2 c 2 , 


By using the knowledge that T - S >0 , the roots of the characteristic equation are obtained from the 
quadratic formula; that is, 




( 14 ) 
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where i=\[-i. Solution of this equation for X yields four roots of equation (13) that are pairs of 
complex conjugates given by 


^"1,2,3 ,4 =± ± / /t+s) 

The homogeneous solution of equation (12) can be written as follows: 

W[[(x) = iqe x sin^T+S x + K 2 ) + K 3 e“ v ^ i= ® (L_X) sin(^T+S x + K 4 ) 


(15) 


(16a) 


where x e [0,L] . The symbols K 1? K 2 , K 3 , and K 4 are real -valued constants that are determined from 
the boundary conditions given by equations (B18) and (B19). The solution given by equation (16a) rep- 
resents a damped, oscillatory response that decays from each end of the cylinder. The regions near the 
edges of the cylinder, where the amplitude of w H (x) is the largest, are called the bending boundary lay- 
ers. All response quantities that exhibit bending boundary layers involve derivatives of equation (16a) 
and can be expressed in the general form 

F(x)=F 1 e _ ^ X sin(/f^x + F 2 )+F 3 e _ ^ (l ^ x) sin( N /TTs x+F 4 ) (16b) 


where F 1 through F 4 are constants. 


When the length of the bending boundary layers is less than half the cylinder length, which is typi- 
cal, equations (16) can be partitioned into one part that applies to the edge x = 0 and the other that 
applies to the edge x = L. The response quantities for the region near x = 0 are obtained by setting 
F 3 = 0 in equation (16b). Similarly, the response quantities for the region near x = L are obtained by 
setting F^ = 0 in equation (16b). 


Formulas for Attenuation Length 

Formulas for the attenuation or decay length of the bending boundary layers are obtained by first 
noting that the response quantities for the region near x = 0 are bounded by the two functions 

±F 1 e _ ^ T_nS X and that the response quantities for the region near x = L are bounded by the two func- 
tions ± F 3 e S ^ x) . Let d denote the length for which the solution attenuates or decays to a value of 
8 times the amplitude F^ or F 3 . A reasonable estimate of the attenuation length or decay length d is 
obtained by replacing x and (L - x) with d in the exponential terms of equation (16b) and by noting 
that the amplitude of F(x) is attenuated by the exponential terms. Thus, the attenuation length or decay 

length is given by e S d = £, which yields 


d =- In e(T - S) -1/2 < — 

2 


( 17 ) 


8 



By using equations (B55) and (11), equation (17) is expressed as 


d d° 

— ,4 

/Rh V'Rh 


( 18 ) 


where d° is the attenuation length, in which anisotropy is neglected {/4— 1), that is, given in nondimen- 
sional form by 


d ° __ ln£ „ 


(19) 


The symbol h is the shell wall thickness, and 0 and /( are nondimensional orthotropy and anisotropy 
parameters or factors, respectively, that are given by 
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(20b) 


where the symbols in these equations are defined in appendix B. 

Other useful fomis of equation (20a) are obtained by introducing an effective membrane Poisson’s 
A12 

ratio v m = — -j= , which is the geometric mean of the two Poisson effects associated with the 

V A 11 A 22 

inplane principal direction of a homogenized ortho tropic material. By using this effective membrane 
Poisson’s ratio, equation (20a) is expressed as 
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( 21 ) 


For a single layer of homogeneous, specially orthotropic material, v m 


7 v 12 v 21 >^ =1 >and 
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E 2 (l -Vi 2 v 2 i) 
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( 22 ) 
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which, when substituted into equation (19), yields results identical to the results presented by Kraus 
(ref. 1), where the decay tolerance is given by 8 = Likewise, for a single layer of isotropic material 
with an arbitrary thickness, v m = v, /i- 1, and 


6 = 


1 

V(1 -v 2 ) 


(23) 


For values of 0.2 < v < 0.33 , 6 varies by less than 2 percent. A 90-percent-decay length (8 = 0. 1) that is 

d° , ^ 

a good approximation to the behavior of homogeneous, metallic shells is given by -- = 1.79 
Applying this formula, for example, to the Space Shuttle solid rocket booster described in references 20 
and 21 (R = 72 in., h = 0.5 in.) gives d° = 0.15R = 10.8 in. 



It is interesting to note that the differences between the attenuation lengths that are based on the 
Sanders-Koiter, the Love-Kirchhoff, and Donne IT s equations appear in the coefficient c and in the 
symbols with overbars in equation (20b) for the anisotropy factor A (see eqs. (B22), (B31)-(B35), and 


(B45)-(B49)). For these equations, the Sanders-Koiter theory is given by |Ll= 


3 

2 


and the Love- 


Kirchhoff theory is given by |A = 1. Donne IF s equations are given by q = 0 . For isotropic and specially 
orthotropic cylinders, /t— 1 and the three sets of shell equations yield identical results. Similarly, for 
antisymmetric cross-ply cylinders (A 16 = A 26 = D 16 = D 26 = B 16 = B 26 = B 12 = B 66 = 0) 
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and the three shell theories yield identical results. 


Simplified Formulas for/* 


For balanced, symmetrically laminated cylinders, the only anisotropic constitutive temis are D 16 


and D 2 6, and the anisotropy factor is given by A = \ / e where 
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For thin-shell theories, such as the Sanders-Koiter theory and the Love-Kirchhoff theory, ^ ^ . This 

result suggests that a useful approximation to equation (25) and the anisotropy factor can be obtained 

from a power series expansion for small values of — . This process yields 
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In this expression, 0<p<— and 0< 


°16 

< 1. Thus, the approximate formula for A indicates 

2 A 66 D n h 2 

that for most practical applications of thin-shell theory, the differences between the three different shell 
theories considered herein and the effect of the flexural anisotropy of a general symmetrically laminated 
cylinder are negligible. 


A simplified formula for the anisotropy factor can be derived for the general expression for A that 
is given by equation (20b). For this case, the following power series expansions for small values of 

R 

are used: 
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Substituting equations (27)-(31) into equation (20b) and expanding the resulting expression in a similar 
manner yields 
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The coefficient is a very complicated expression, and as a result, the following first-order approxi- 
mation of /4 is used herein; that is, 


*4 1 + 



(33) 


where /i ^ is the value of equation (20b) with p = 0, which is the anisotropy factor that corresponds to 
the use of Donnell’s equations. The expression for is given by 



In this expression, a 22 and b 21 are obtained from equations (B45)-(B48) by setting p = 0 in equa- 
tions (B31)-(B35). The expression for is obtained from equation (B49) in a similar manner. The 
temi Si in equation (33) represents a first-order correction to the results that correspond to DonnelTs 
equations and is given by 


a 22 D ll e 0 a 22 e l +a 122^0 


b 121 e 0~ b 21 e l -2a 122 b 21 e 0 


4a 22^0 \/ a 22 D ll e 0 -b 21 1 


where the temis that appear in equation (35) are given in appendix C. In addition, further simplifications 
to an d C\ are a l so presented in appendix C for unbalanced and balanced symmetric laminates and 
for balanced, unsymmetric laminates that include the subclasses of general antisymmetric laminates, 
antisymmetric cross-ply laminates, and antisymmetric angle-ply laminates. The relative size of S\ and 
its contribution to equation (33) are examined parametrically in the subsequent section of the present 
study. 

Results and Discussion 

Equations (18) and (19) form the basis for the parametric study presented herein. In particular, 
the two equations isolate the contributions of shell orthotropy and shell anisotropy to the bending 
boundary-layer decay length with nondimensional parameters and imply the generic design-chart repre- 
sentations that are illustrated in figures 1 and 2. In figure 1, generic results are presented that show the 

nondimensional, 90-percent-decay length given by ^ , as a function of the orthotropy param- 

£ = 0.1 

eter 0 , for selected values of the anisotropy parameter /t. A 90-percent-decay length was selected herein 
to yield an accuracy that is approximately to within the accuracy of the experimentally detemiined 
material properties, but other values could be used. 

In a manner similar to figure 1, figure 2 shows the nondimensional, 90-percent-decay length, as a 
function of the anisotropy parameter /t, for selected values of the orthotropy parameter 0. Results that 
correspond to balanced, symmetrically laminated cylinders are given by a value of /l- 1, and results 
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that correspond to an isotropic shell wall are indicated in the figures by the filled circle with an ordinate 
value of 1.79. Overall, these two figures represent results that are applicable to a vast range of laminate 
constructions and provide a common basis for comparison of regular and hybrid laminates made of dif- 
ferent material systems and laminate stacking sequences. In general, the figures show increases in the 
nondimensional 90-percent-decay length, with increases in either the orthotropy parameter 0 or the 
anisotropy parameter A. In addition, the results in figures 1 and 2 clearly indicate the effect of neglect- 
ing shell-wall anisotropy on the attenuation of a bending boundary layer. 


The actual value of the nondimensional, 90-percent-decay length depends on the particular values 
of the orthotropy and anisotropy parameters of a given laminate. Thus, additional results are presented 
subsequently that show how the orthotropy parameter 0 and the anisotropy parameter /t vary with 
laminate construction. In particular, values of 0 and /4 are presented first for balanced and unbalanced 
symmetrically laminated cylinders. Then, values are presented for balanced and unbalanced unsymmet- 
rically laminated cylinders. Nine different contemporary material systems were used to generate 
these results. These material systems include boron-aluminum, S-glass-epoxy, a typical boron-epoxy, 
AS4/3501-6 graphite-epoxy, AS4/3502 graphite-epoxy, IM7/5260 graphite-bismaleimide, Kevlar 
49-epoxy, IM7/PETI-5 graphite-epoxy, and P-100/3502 pitch-epoxy materials. The mechanical 
properties of these material systems are presented in table 1; the nominal ply thickness that was used is 
0.005 in. An example problem in appendix D demonstrates how to perfomi the calculations that were 
used to generate the results that are presented subsequently. 


Balanced, Symmetrically Laminated Cylinders 

Symmetrically laminated shell walls are characterized mathematically by values of zero for the sub- 
scripted B temis that appear in the constitutive equation (A 15). In addition, balanced, symmetrically 
laminated shell walls do not exhibit coupling between extension and shear, which is characterized by 
Ai6 = A 26 = 0 in equation (A15). Shell walls of this class are strictly specially orthotropic for many 
laminates. However, for some wall constructions, balanced, symmetric laminates exhibit anisotropy in 
the fomi of coupling between pure bending and twisting of the shell wall. This type of anisotropy is 
manifested by nonzero values of the D 16 and D 26 constitutive terms in equation (A15). However, the 
discussion of equation (26) that has been given herein indicates that this type of anisotropy is negligible 
for thin shells and that the differences between results obtained from the Sanders-Koiter, the 
Love-Kirchhoff, and the Donnell theories are insignificant. Moreover, — 1 for this class of laminated- 
composite shell walls, and the attenuation behavior is governed by the nondimensional orthotropy 
parameter 0 that is given by equation (20a). Furthermore, equations (18) and (19) indicate that the 
attenuation length is a constant multiple of the orthotropy parameter that depends on the attenuation- 
tolerance parameter e. For this case, trends that are exhibited by 0 are identical to those exhibited by the 
attenuation length based on any value of 8. 


Values of the orthotropy parameter 0 are presented in figure 3 and table 2 for single-ply, homoge- 
neous, specially orthotropic and isotropic shell walls, with arbitrary thickness, as a function of the 
ratio of the principal elastic moduli, E 2 /Ei. For these results, the orthotropy parameter is given by 
equation (22) and is expressed in the following, more convenient form: 


0 =\ 
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One curve (shown in fig. 3), which is essentially several coincident curves, corresponds to general 
results for 0.2 <v 12-° .35. In addition, specific results for the nine material systems considered herein 
and for a typical aluminum and a steel are indicated by the square symbols in the figure. The results in 
figure 3 indicate that the effects of variations in the major Poisson’s ratio on the ortho tropy parameter 0 
are small compared to the effect of variations in the ratio of the principal elastic moduli. Moreover, the 
results show that 0 decreases rapidly as the ratio of the principal elastic moduli increases, particularly 
for values of E 2 /Ei less than approximately 0.1, which corresponds to most of the contemporary ortho- 
tropic materials considered herein. Figure 3 also shows that an isotropic material corresponds to 0 ~ 1. 

Values of the orthotropy parameter 0 for the single-ply, homogeneous, specially orthotropic cylin- 
ders investigated by Cheng and He (ref. 6) were also obtained. A comparison of the results obtained in 
the present study, with the corresponding results of reference 6, is presented in table 3 for boron-epoxy, 
glass-epoxy, and graphite-epoxy materials and for the cylinder radius-to-thickness ratio R/h = 208.311. 
Moreover, a range of results is shown for reference 6 that corresponds to various simplifications that 
were used in the equations that govern the response. The actual material properties that were used are 
given in reference 6. In table 3, the quantity used for comparison is given by 



which is the real part of the exponent that appears in the eigenfunction solution used by Cheng and 
He (n = 0 in eq. (25) of ref. 6; see also eq. (47) of ref. 13), which corresponds to the decay or attenuation 
of the response. The orthotropy parameter shown in equation (37) is defined by equation (36). The 
results in table 3 show very good agreement (less than 1 percent difference) for all three materials. In 
addition, the results obtained herein that are shown in table 3 for the boron-epoxy material are also in 
excellent agreement with the corresponding results presented by McDevitt and Simmonds (ref. 13). 

Values of the orthotropy parameter 0 are presented in figure 4 and table 4 for multilayered [(±(|)) m ] s 
laminates made from the nine material systems as a function of the fiber angle 0, which is measured 
from the x-axis toward the 0-axis. The results are independent of the stacking sequence number m and 
show a wide variation in 0 with the material system. The results also show, for the most part, a wide 
variation in 0 with the fiber angle § and a reduction in 0 as the fiber angle increases from 0° to 90°. The 
largest value (2.93) and the smallest value (0.34) of 0 are exhibited by the unidirectional laminates 

made from P-100/3502 pitch-epoxy material and correspond to values of ^ equal to 5.13 and 

£ = o.l 

0.59, respectively. Moreover, the greatest variation in 0 with the fiber angle (approximately 8.7 times) 
is exhibited by the laminates made from P-100/3502 pitch-epoxy material. The smallest variation is 
exhibited by the laminates made from the boron-aluminum material. 

Results are presented in figure 5 that show the values of the orthotropy parameter for [(±45/0 2 ) m ] s , 
[(0 2 /±45) m ] s , [(±45/90 2 )J s , [(90 2 /±45)J s , [(±45/0/90) m ] s , and [(0/90/±45) m ] s laminates made of 
IM7/5260 graphite-bismaleimide material for values of the stacking sequence number m = 1 to 6. 
Values of 0 range from approximately 1.53 to 0.64. These results show that the curves for the 
[(±45/0 2 ) m ] s and [(0 2 /±45) m ] s laminates approach 0 ~ 1.41 as m increases to a value of 6, with the 
curve for the [(0 2 /±45) m ] s laminates converging from above and the other curve converging from 
below. The higher values of 0 for the [(0 2 /±45) m ] s laminates are attributed to the higher axial bending 
stiffness that is obtained by placing the 0° plies at the outer surfaces of the laminates, particularly for the 
lower values of the stacking sequence number m. Similarly, the results in figure 5 show that the curves 
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for the [(±45/90 2 ) m ] s and [(902/±45) m ] s laminates approach 0 -0.76 as m increases to a value of 6, 
with the curve for the [(±45/902) m ] s laminates converging from above and the other curve converging 
from below. Likewise, the results in figure 5 show that the curves for the [(±45/0/90) m ] s and 
[(0/90/±45) m ] s quasi-isotropic laminates approach 0 - 1.03 as m increases to a value of 6, with the 
curve for the [(0/90/±45) m ] s laminates converging from above and the other curve converging from 
below. 

Overall, the results in figure 5 indicate that the [(±45/0 2 ) m ] s and [(0 2 /±45) m ] s laminates exhibit 
higher values of the ortho tropy parameter than the [(±45/0/90) m ] s and [(0/90/±45) m ] s quasi-isotropic 
laminates, which exhibit higher values of the orthotropy parameter than the [(±45/90 2 ) m ] s and 
[(90 2 /±45 ) m ] s laminates. This trend corresponds to a reduction in the value of 0 as the axial bending 
and extensional stiffnesses of the laminates decrease. 

Results are presented in figure 6 and table 5 that show the effect of the nine material systems con- 
sidered herein on the orthotropy parameter for the [(0 2 /±45) m ] s laminates. Values of 0 range from 
approximately 1.67 for P-100/3502 pitch-epoxy material to 1.09 for boron-aluminum material. Most of 
the materials exhibit values of 6 in the range of approximately 1.4 to 1.6. All curves show about the 
same reduction in 6 as the stacking sequence number m increases. 

Results similar to those in figure 6 and table 5 are presented in figure 7 and tables 6 and 7 that 
show the effect of the nine material systems on the orthotropy parameter for the [ (±45/0/90) m ] s and 
[(0/90/±45) m ] s quasi-isotropic laminates. These results show a much smaller variation in the orthotropy 
parameter with material system and stacking sequence number for the quasi-isotropic laminates than for 
the [(02/±45) m ] s laminates in figure 6. In particular, values of 6 for the quasi-isotropic laminates range 
from approximately 1. 15 to 1.0. The largest values of 6 in figure 7 are exhibited by the laminates made 
of the P-100/3502 pitch-epoxy material. Moreover, the results show a larger variation in 0 with stack- 
ing sequence number for the [(0/90/±45) m ] s laminates than for the [(±45/0/90) m ] s laminates. 

Unbalanced, Symmetrically Laminated Cylinders 

Unbalanced, symmetric laminates exhibit anisotropy in the fomi of extensional-shear coupling 
(A 16 ^ 0, A 2 6 ^ 0) in addition to flexural anisotropy (D 16 ^ 0, D 2 £ ^ 0). For these laminates, the value of 
the anisotropy parameter , given by equations (20b) and (33), is not equal to unity. Simplified 
expressions for the anisotropy parameter and the first-order correction factor defined by equa- 
tions (33)-(35), are given by equations (C23) and (C24), respectively. Equation (C24) indicates that the 
value depends on coupling between the membrane and flexural anisotropies. 

Values of the orthotropy parameter 0 for [(+(|)) 2 m ]s symmetric, unidirectional laminates for the 
nine material systems considered herein are also presented in figure 4 and table 4; that is, the curves 
presented in figure 4 and the data presented in table 4 for the [(±4>) m ] s symmetric angle-ply laminates 
are identical to those for the corresponding [(+402nJs symmetric, unidirectional laminates. Thus, the 
orthotropy behavioral trends for the undirectional laminates are identical to those discussed previously 
for the symmetric angle-ply laminates and are also independent of the stacking sequence number m. 

Results for the anisotropy parameter /t ^ and the first-order correction factor are shown in fig- 
ure 8 and table 8 and in figure 9, respectively, for the IX+^nJs symmetric, unidirectional laminates 
with the nine material systems considered herein and are independent of the stacking sequence number 
m. The results in figure 8 and table 8 show a substantial variation in ^ with fiber orientation and with 
material system. The results show that A ^ is the most pronounced for values of the fiber angle § 
between approximately 55° and 80° and that the contribution of the anisotropy to the attenuation 
behavior is essentially insignificant (<1.05) for values of § < 25° and 0 > 85°. Moreover, the largest 
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variation in At Q with fiber angle is exhibited by the laminates made of the P-100/3502 pitch-epoxy 
material, and the smallest variation is exhibited by the laminates made of boron-aluminum material. 
Values of A range from approximately 1.42 for the maximum point on the curve for the P-100/3502 
pitch-epoxy material to a value of 1 . 


The results shown in figure 9 for the first-order correction factor Si for the [(+(|))2 m ]s symmetric, 
unidirectional laminates indicate a substantial relative variation in S ^ with fiber orientation and with 
material system, but all values of Si are less than approximately 0.45. Moreover, Si is less than 
approximately 0.2 for all materials except the P-100/3502 pitch-epoxy material. For the upper bound of 

thinness of thin-shell theory, given by ^ the contribution of S\ to the anisotropy factor 

defined by equation (33) is practically negligible. Equation (C24) indicates that the insignificance of 
Si means that the coupling of the membrane and flexural anisotropies is negligible for these laminates. 


The insignificance of Si is illustrated and verified in figure 10 for the [(+(|)) 2m ] s symmetric, unidirec- 
tional laminates made of IM7/5260 graphite-bismaleimide material (black curves) and P-100/3502 

pitch-epoxy material (gray curves) for ^ The finely dashed curves shown in figure 10 

R 20 

correspond to 90-percent-decay lengths for which the anisotropy is neglected. In contrast, the solid 
curves and the coarsely dashed gray curve include the effect of the membrane anisotropy and are shown 
for values of g = 0, 1.0, and 1.5. For these values, results that correspond to the Sanders-Koiter theory 

3 

and the Love-Kirchhoff theory are given by g = - and g = 1.0, respectively. Results that correspond to 

Donnell’s equations are given by g = 0. The solid curves in figure 10 for g . = 1.0 and 1.5 are based on 
the exact solution that uses equation (20b) for the anisotropy factor. The corresponding curves that are 
based on the approximate formula for the anisotropy parameter that is given by equation (33) are identi- 
cal. The solid curves and the coarsely dashed gray curve indicate that varying g yields a small effect, 
which implies that all three shell theories yield essentially the same results and that At ~ At ^ for the 
[(+(|)) 2m ] s symmetric, unidirectional laminates. Comparing the solid and finely dashed curves in 
figure 10 also indicates that neglecting the membrane anisotropy underestimates the bending boundary- 
layer decay length by as much as approximately 31 percent and 21 percent for shell walls made of 
P-100/3502 pitch-epoxy and IM7/5260 graphite-bismaleimide materials, respectively. 


Values of the ortho tropy parameter 0 for [(+45 2 /0/90) m ] s and [(0/90/+45 2 ) m ] s laminates made of 
the nine material systems considered herein are also presented in figure 7 and in tables 6 and 7, respec- 
tively. More specifically, the values of 0 for these laminates are identical to the values for the corre- 
sponding quasi-isotropic laminates. Results for the anisotropy parameter defined by equation (34), 
are shown in figure 11 for [(+45 2 /0/90) m ] s and [(0/90/-h45 2 ) m ] s laminates made of the nine material 
systems considered herein. The results in figure 11 are identical for the two laminate families, show no 
significant variation in At ^ with the stacking sequence number m, and only a slight variation (less than 
approximately 9 percent) with the material system. Values of At§ range between approximately 1.1 
and 1.0. Corresponding results for the first-order correction factor S^ defined by equation (35), 
which are not shown herein, were obtained and indicate that all values of Si for the [(+45 2 /0/90) m ] s 
and [(0/90/+45 2 ) m ] s laminates are less than approximately 0.1. These values indicate that the contribu- 
tion of Si to the anisotropy factor defined by equation (33) is practically negligible. Thus, At ~At^ for 
these laminates. The values of /t ^ shown in figure 11 suggest that neglecting the anisotropy would, at 
most, underestimate the bending boundary-layer decay length by approximately 10 percent. The insig- 
nificance of Si also means that the coupling of the membrane and flexural anisotropies are unimportant 
with regard to the primary effect of the individual shell anisotropies that is captured by the parameter 
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Balanced, Unsymmetrically Laminated Cylinders 

Balanced, unsymmetric laminates may, in general, exhibit anisotropy in the form of coupling 
between pure bending and twisting (D 16 =£ 0, D 26 0) and coupling between membrane and bending 

action, which is manifested by nonzero values for any of the subscripted B-terms in equation (A15). 
These laminates do not, however, exhibit extensional-shear coupling (A^ = A 2 £ = 0). For the unsym- 
metric laminates that are discussed subsequently, the first ply in the stacking sequence is the innermost 
ply of a cylinder. Simplified expressions for the anisotropy parameter Aq and the first-order correction 
factor defined by equations (33)-(35), are given by equations (C27) and (C29), respectively. Equa- 
tions (C28) and (C29) indicate that the value of ^ depends on coupling between the flexural anisot- 
ropy and the anisotropy caused by unsymmetric lamination. 

Results for regular, antisymmetric angle-ply laminates are shown in figures 4 and 12-16, and in 
tables 4 and 9. In particular, values of the ortho tropy parameter 0 for [(d=c])) m ] T unsymmetric laminates 
made of the nine material systems considered herein are also presented in figure 4 and table 4; that is, 
the ortho tropy-parameter curves presented in figure 4 and the data presented in table 4 for the [(±(|)) m ] s 
symmetric angle-ply laminates are also identical to those for [(=tc[)) m ] T unsymmetric laminates. Thus, the 
orthotropy behavioral trends for the [(=tc|)) m ] T unsymmetric laminates are identical to those discussed 
previously for the corresponding symmetric angle-ply laminates and are also independent of the stack- 
ing sequence number m. 

Results for the anisotropy parameter A ^ defined by equation (C27) are shown in figure 12 and 
table 9 for two-ply [±cj)] T unsymmetric laminates made of the nine material systems considered herein. 
The results in figure 12 show a substantial variation in Aq with fiber orientation and with material 
system and also show that A q is the most pronounced for values of the fiber angle 0 between approxi- 
mately 15° and 60°. Moreover, the largest variation in Aq with fiber angle is exhibited by the laminates 
made of the P-100/3502 pitch-epoxy material, and the smallest variation is exhibited by the laminates 
made of boron- aluminum material. Values of Aq range from approximately 0.75 for the minimum 
point on the curve for the P-100/3502 pitch-epoxy material to a value of 1.0. The results in figure 13 
show the variation in A 0 with the fiber angle § and the stacking sequence number m for | (±(|)) m ] T 
unsymmetric laminates made of the P-100/3502 pitch-epoxy material. These results show a rapid 
decline in the importance of Aq that is manifested by the curve moving closer to Aq = 1, as the stack- 
ing sequence number increases. For m = 2, 0.95 < Aq < 1. 

Results for the first-order correction factor are shown in figure 14 for two-ply [±(|)] T unsymmet- 
ric laminates made of the nine material systems considered herein. The results in figure 14 also show a 
substantial variation in with fiber orientation and with material system. However, the maximum 
value of dq < 0.07 for all material systems. Results are presented in figure 15 that show the variation in 
&Y with the fiber angle 0 and the stacking sequence number m for [ (±(f)) m ] T unsymmetric laminates 
made of the P-100/3502 pitch-epoxy material. These results show significant reductions in with an 
increase in the stacking sequence number. 

Overall, the results in figures 14 and 15 indicate that the contribution of to the anisotropy factor 
defined by equation (33) is negligible for the upper bound of thinness given by — = — , which means 
that A ~Aq. Thus, the results in figure 12 for the two-ply [±(|)] T unsymmetric laminates indicate that 
neglecting the shell anisotropy overestimates the bending boundary-layer decay length (because 
A 0 < 1) by as much as approximately 33 and 22 percent for shell walls made of P-100/3502 pitch- 
epoxy and IM7/5260 graphite-bismaleimide materials, respectively. The insignificance of ^ is illus- 
trated in figure 16 by the gray and the black curves for the laminates made of P-100/3502 pitch-epoxy 
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and IM7/5260 graphite-bismaleimide materials, respectively. The solid black and the gray curves are 
for the upper bound of thin-shell theory that is given by — = The finely dashed curves shown in 
figure 16 correspond to 90-percent-decay lengths for which the anisotropy is neglected. In contrast, the 
solid curves include the effect of the shell anisotropy and are shown for values of p = 0, 1.0, and 1.5. 
The solid curves for p . = 1.0 and 1.5 are based on the exact solution that uses equation (20b). The corre- 
sponding curves that are based on the approximate formula for the anisotropy parameter given by 
equation (33) are identical. The solid curves indicate no significant effect of varying p, which implies 
that all three shell theories yield essentially the same results for the [±c])] T unsymmetric laminates. For 
[(±(|)) m ] T unsymmetric laminates with m > 1 that are made from any of the nine material systems 
considered herein, the results in figures 12 through 15 indicate that neglecting the shell- wall anisotropy 
will have a small effect on the calculation of the bending boundary-layer decay length. 

Values of the orthotropy parameter 6 and the anisotropy parameter for [0p/90 q ]x unsymmetric 
cross-ply laminates are shown in figure 17 and table 10, and in figure 18 and table 11, respectively, for 
the nine material systems considered herein and as a function of the percentage of 0° plies. For this class 
of laminates, equation (20b) simplifies to equation (34); that is, /t =s?q. This simplification means that 
the anisotropy parameter is independent of p, which means that all three shell theories considered 
herein yield identical results. 

The results in figure 17 show a large variation in 0 with the percentage of 0° plies for most of the 
material systems. In addition, the results show a large variation in 6 with the material system for the 
laminates that are dominated by 90° plies (less than approximately 10 percent 0° plies) and by 0° plies 
(more than approximately 80 percent 0° plies). Values of 0 vary the most for laminates made of 
P-100/3502 pitch-epoxy material, with values that range from approximately 0.3 to 2.93. Most of the 
materials exhibit values of 0 in the range of approximately 0.6 to 2.1. 

The results in figure 18 also show a large variation in with the percentage of 0° plies for most 
of the material systems and a large variation with material system for laminates with <70 percent 
0° plies. Moreover, the results show that A ^ is the most pronounced (most different from a value of 1) 
for laminates with approximately 15 to 30 percent 0° plies. The largest variation in A ^ with the per- 
centage of 0° plies is exhibited by the laminates made of the P-100/3502 pitch-epoxy material, and the 
smallest variation is exhibited by the laminates made of boron-aluminum material. Values of range 
from approximately 0.57 for the minimum point on the curve for the P-100/3502 pitch-epoxy material 
to a value of 1. Thus, in some cases neglecting the shell-wall anisotropy overestimates the bending 
boundary-layer decay length by as much as approximately 75 percent for a shell wall made of 
P-100/3502 pitch-epoxy material. This result is illustrated in figure 19 by the gray curves. Similar 
results are presented in figure 19 for [0p/90 q ] T unsymmetric cross-ply laminates made of IM7/5260 
graphite-bismaleimide material (black curves). The solid black and the gray curves include the effect of 
the shell anisotropy, and the finely dashed curves shown in the figure correspond to 90-percent-decay 
lengths for which the anisotropy is neglected. The results in figure 19 show that including the effect of 
anisotropy is particularly important for laminates with less than approximately 70 percent 0° plies and 
more than approximately 5 percent 0° plies. 

Unbalanced, Unsymmetrically Laminated Cylinders 

Unbalanced, unsymmetric laminates may, in general, exhibit full anisotropy in the fomi of coupling 
between pure bending and twisting (D 16 ^ 0, D 2 6 ^ 0) and coupling between membrane and bending 
action, which is manifested by nonzero values for any of the subscripted B -terms in equation (A15), and 
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extensional-shear coupling (A 16 ^ 0, A 2 6 ^ 0). The expressions for the anisotropy parameter A^ and 
the first-order correction factor that are given by equations (C2)-(C22) indicate that A ^ exhibits 
coupling between the membrane anisotropy and the anisotropy that is caused by unsymmetric lamina- 
tion, and that exhibits coupling between all three types of anisotropies. One family of laminates that 
exhibits all these anisotropies is the [70p/0 q ] T unbalanced, unsymmetric laminates with p > 0 and q ^ 0. 

Values of the orthotropy parameter 0 and the anisotropy parameter A ^ for [70p/0 q ]x unbalanced, 
unsymmetric laminates are shown in figure 20 and table 12 and in figure 21 and table 13, respectively, 
for the nine material systems considered herein and as a function of the percentage of 70° plies. The 
results in figure 20 show a large variation in 0 , with the percentage of 70° plies for most of the material 
systems. The results also show a large variation in 0 with the material system for the laminates that 
aredominated by 0° plies (less than approximately 20 percent 70° plies). Values of 0 vary the most for 
laminates made of P-100/3502 pitch-epoxy material, with values that range from approximately 0.5 
to 3.0. 

The results in figure 21 also show a substantial variation in A ? 0 with the percentage of 70° plies for 
most of the material systems and a large variation with material system for laminates with between 
approximately 45 and 100 percent 70° plies. The largest overall variation in A ^ with the percentage of 
70° plies is exhibited by the laminates made of the P-100/3502 pitch-epoxy material, and the smallest 
variation is exhibited by the laminates made of boron- aluminum material. Values of A ^ range from 
approximately 1.4 to 0.95, which correspond to the maximum and minimum points, respectively, on the 
curve for the P-100/3502 pitch-epoxy material. 


Results for the first-order correction factor were also obtained for [70 p /0 q ]x unbalanced, 
unsymmetric laminates made of the nine material systems considered herein but are not included in the 
present paper. These results also show a substantial, relative variation in ^ with the percentage of 

70° plies, but overall, the magnitude of is less than approximately 0.25 for the P-100/3502 pitch- 

epoxy material and <0.1 for the other materials. These results indicate that the contribution of ^ to the 
anisotropy factor defined by equation (33) is negligible for the upper bound of thin-shell theory that is 

given by ^ which means that A~A q. Thus, the results in figure 21 suggest that, in some 

cases, neglecting the shell- wall anisotropy may overestimate the bending boundary-layer decay length 
and, in other cases, may underestimate the decay length. This statement is based on the observation that, 
for example, when Aq < 1 , including anisotropy reduces the value of the decay length given by 

equations (18) and (19). The insignificance of also means that the contribution of the flexural anisot- 
ropy to the coupling of the anisotropies is negligible. The insignificance of is clarified in figure 22 

for laminates made of P-100/3502 pitch-epoxy material (gray curves) and of IM7/5260 graphite- 
bismaleimide material (black curves). The solid black and the gray curves are for the upper bound of 

thinness given by — = The finely dashed curves shown in the figure correspond to 90-percent- 
R 20 

decay lengths for which the anisotropy is neglected. In contrast, the solid curves include the effect of the 
shell anisotropy and are shown for values of p = 0, 1.0, and 1.5. Moreover, the solid curves for p , = 1.0 
and 1.5 are based on the exact solution that uses equation (20b). The corresponding curves that are 
based on the approximate formula for the anisotropy parameter that is given by equation (33) are identi- 
cal. The solid curves indicate a negligible effect of varying p, which verifies that A ~A 0 and implies 
that all three shell theories yield essentially the same results for the [70 p /0 q ] T unbalanced, unsymmetric 
laminates. In addition, the results show that neglecting the shell-wall anisotropy, for the most part, 
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underestimates the bending boundary-layer decay length by as much as approximately 16 and 6 percent 
for shell walls made of P-100/3502 pitch-epoxy and IM7/5260 graphite-bismaleimide materials, 
respectively, and with approximately 20-percent 70° plies. In addition, the results in figure 22 show that 
neglecting the shell-wall anisotropy underestimates the bending boundary-layer decay length by as 
much as approximately 31 and 20 percent for shell walls made of P-100/3502 pitch-epoxy and 
IM7/5260 graphite-bismaleimide materials, respectively, and with approximately 100-percent 70° plies. 
There is only a very small range shown in figure 22 where neglecting the shell-wall anisotropy overesti- 
mates the bending boundary-layer decay length, and for this region, the effect is negligible. 


Concluding Remarks 

An analytical study of the attenuation of bending boundary layers in both balanced and unbalanced, 
symmetrically and unsymmetrically laminated-composite, thin cylindrical shells has been presented for 
nine contemporary material systems. The analysis is based on the linear Sanders-Koiter shell equations 
and contains the Love-Kirchhoff shell equations and DonnelTs equations as special cases. With this 
analysis, two nondimensional parameters have been identified that characterize and quantify the effects 
of laminate orthotropy and laminate anisotropy on the bending boundary-layer decay length in a very 
general and encompassing manner. The anisotropy parameter includes the effects of anisotropy in the 
fomi of coupling between pure bending and twisting that appears in many symmetric laminates to some 
extent, coupling between extension and shear that is present in unbalanced laminates, and coupling 
between membrane and bending action that is present in unsymmetric laminates. 

A substantial number of structural design technology results for the bending boundary-layer decay 
length have been presented for a wide range of laminated-composite shell structures that should be use- 
ful additions to the structural designer’s collection of preliminary design tools. Moreover, the analysis 
and results should provide additional physical insight into the fundamental behavior of general 
laminated-composite shell structures and create a common basis for assessing bending boundary-layer 
attenuation for the vast range of laminate constructions that are possible. Furthermore, the results 
should be useful for the design of specimens for material characterization tests, for instrumenting struc- 
tural verification tests, and for defining finite-element meshes. For all the laminate constructions con- 
sidered in the present study, the results show that the differences between results that were obtained 
with the Sanders-Koiter shell equations, the Love-Kirchhoff shell equations, and DonnelTs equations 
are negligible. The results also show that the effect of anisotropy in the form of coupling between pure 
bending and twisting has a negligible effect on the size of the bending boundary-layer decay length of 
the balanced, symmetrically laminated cylinders considered. Moreover, the results show that the cou- 
pling of the membrane and flexural anisotropy and the anisotropy caused by unsymmetric lamination is 
generally unimportant with regard to the primary effect of the individual shell anisotropies on the bend- 
ing boundary-layer decay length. The only exception encountered was for unbalanced, unsymmetrically 
laminated cylinders for which coupling of the membrane anisotropy and the anisotropy caused by 
unsymmetric lamination is a primary effect, as expected. The results also show that in some cases 
neglecting the shell anisotropy results in underestimating the bending boundary-layer decay length and 
in other cases results in an overestimation. 
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Table 1. Lamina Properties 


Lamina 

property 

(a) 

Material systems 

Boron- 

A1 

S -glass- 
epoxy 

Kevlar 

49-epoxy 

IM7/ 

5260 

AS4/ 

3502 

AS4/ 

3501-6 

Boron- 

epoxy 

IM7/ 

PETI-5 

P-100/ 

3502 

E 1? Msi 

33 

7.5 

11.02 

22.1 


20.01 

29.58 

20.35 

53.5 

E 2 , Msi 

21 

1.7 

0.8 

1.457 


1.30 

2.68 

1.16 

0.73 

v 12 

0.23 

0.25 

0.34 

0.258 


0.30 

0.23 

0.29 

0.31 

Gj 2 , Msi 

7.0 

0.80 

0.33 

0.860 

0.87 

1.03 

0.81 

0.61 

0.76 

OCJ x 10 6 /°F 

3.2 

3.5 

-2.22 

0.0125 

0.25 

-0.167 

3.38 

-0.14 

-0.64 

a 2 x 10 6 /°F 

11.0 

11.0 

43.89 

14.91 

16.2 

15.6 

16.83 

16.85 

17.2 


^Subscripts 1 and 2 denote the longitudinal (fiber) and transverse (matrix) directions of a specially orthotropic 
lamina, respectively. 


Table 2. Orthotropy Parameters for Single-Ply, Homogeneous, 
Specially Orthotropic Laminates 

[See eq. (36)] 


E 2 /Ei 

Orthotropy parameter, 0 

V 12 — 0.2 

V 12 = 0.35 

0.01 

3.16 

3.16 

0.02 

2.66 

2.66 

0.04 

2.24 

2.24 

0.06 

2.02 

2.02 

0.08 

1.88 

1.88 

0.1 

1.78 

1.78 

0.15 

1.61 

1.61 

0.2 

1.50 

1.50 

0.3 

1.36 

1.36 

0.4 

1.26 

1.27 

0.5 

1.20 

1.21 

0.6 

1.14 

1.16 

0.7 

1.10 

1.12 

0.8 

1.07 

1.09 

0.9 

1.04 

1.06 

1.0 

1.01 

1.03 


Table 3. Results for Specially Orthotropic Materials 


Material 

systems 

E 2 /E! 

(a) 

v 12 

6 

Ke(fi), from 
reference 6 

(b) 

from 
present study 

(b) 

Boron-epoxy 

wmsssmm 

0.30 

1.782 

2.796-2.805 

2.806 

Glass-epoxy 


0.25 

1.323 

3.757-3.779 

3.779 

Graphite -epoxy 

— 

0.25 

2.516 

1.984—1.987 

1.987 


^Subscripts 1 and 2 denote the major and minor principal directions, respectively, of the specially orthotropic 
materials defined in reference 6. 

^Quantity 'ge(fi) is defined by equation (37). 
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Table 4. Orthotropy Parameters for [(±4>) m ] s , and [ ( Laminates 



Orthotropy parameter, 0 

IM7/ 

5260 

Boron- 

A1 

S -glass- 
epoxy 

Kevlar 
49 -epoxy 

AS4/ 

3502 

AS4/ 

3501-6 

Boron- 

epoxy 

IM7/ 

PETI-5 



1.98 

1.13 

WEM 

1.93 

mesm 

1.98 

1.82 

2.05 

2.93 

5 

1.97 

1.13 

mm 

1.92 

litl 

1.97 

1.82 

2.04 

2.90 


1.94 

1.13 

mm 

1.91 

wEm 

1.94 

1.81 

2.02 

2.84 

15 

1.90 

1.12 

mm. 

1.88 

1.77 

1.88 

1.79 

1.98 

2.73 


1.84 

1.12 

mm 

1.84 

1.72 

1.80 

1.76 

1.92 

2.58 

25 

1.76 

1.11 

WSm 

1.78 

1.66 

1.71 

1.72 

1.84 

2.40 

30 

1.65 

1.10 

1.30 

1.70 

1.57 

1.59 

1.67 

1.74 

2.19 

35 

1.53 

1.09 

1.24 

1.60 

1.46 

1.46 

1.59 

1.61 

1.97 

40 

1.39 

1.07 

1.16 

1.46 

1.34 

1.32 

1.48 

1.46 

1.74 

45 

1.23 

1.05 

1.08 

1.30 

1.20 

1.17 

1.33 

1.29 

1.50 

50 

1.07 

1.02 

1.00 

1.13 

1.05 

1.03 

1.15 

1.12 

1.26 

55 

0.92 

0.99 

0.92 

0.96 

0.92 

0.89 

0.97 

0.94 

1.04 

60 

0.78 

0.97 

0.85 

0.81 

0.80 

0.77 

0.81 

0.79 

0.83 

65 

0.68 

0.95 

0.79 

0.69 

0.70 

0.68 

0.70 

0.67 

0.65 

70 

0.60 

0.93 

0.75 

0.61 

0.63 

0.60 

0.62 

0.59 

0.51 

75 

0.55 

0.92 

0.72 

0.56 

0.59 

0.55 

0.58 

0.53 

0.42 

80 

0.52 

0.91 

0.71 

0.53 

0.56 

0.52 

0.56 

0.50 

0.37 

85 

0.51 

0.90 

0.70 

0.52 

0.55 

0.51 

0.55 

0.49 

0.35 

90 

0.51 

0.90 

0.69 

0.52 

0.55 

0.51 

0.55 

0.49 

0.34 


Table 5. Orthotropy Parameters for [(0 2 /±45) m ] s Laminates 


Number of 

Orthotropy parameter, 6 

laminate 

IM7/ 

Boron- 

S -glass- 

Kevlar 

AS4/ 

AS4/ 

Boron- 

IM7/ 


plies, 8m 

5260 

A1 

epoxy 

49-epoxy 

3502 

3501-6 

epoxy 

PETI-5 


8 

1.53 

1.12 

1.32 

1.54 

1.49 

1.52 

1.52 

1.56 

1.67 

16 

1.48 

1.10 

1.28 

1.48 

1.44 

1.46 

1.46 

1.50 

1.60 

24 

1.46 

1.10 

1.27 

1.46 

1.42 

1.44 

1.44 

1.48 

1.58 

32 

1.45 

1.09 

1.26 

1.45 

1.41 

1.43 

1.43 

1.47 

1.56 

40 

1.44 

1.09 

1.26 

1.44 

1.40 

1.42 

1.43 

1.46 

1.56 

48 

1.43 

1.09 

1.25 

1.44 


1.42 

1.42 

1.46 

1.55 

56 

1.43 

1.09 

1.25 

1.44 


1.42 

1.42 

1.45 

1.55 

64 

1.43 

1.09 

1.25 

1.43 


1.41 

1.42 

1.45 

1.54 


23 



































Table 6. Orthotropy Parameters for [(±45/0/90 ) in ] s and [(+45 2 /0/90) m ] s Laminates 


Number of 

Orthotropy parameter, 0 

laminate 
plies, 8m 

IM7/ 

5260 

Boron- 

A1 

S -glass- 
epoxy 

Kevlar 

49-epoxy 

AS4/ 

3502 

AS4/ 

3501-6 

Boron- 

epoxy 

IM7/ 

PETI-5 


8 

LOO 

1.01 

1.01 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

16 

mm, 

mm 

mm 

mm 

mm 

■EK 

mSm 

mm 

mm 

24 

HS 

H 

■H 

ip 

WBm. 

mm 

mm 

m m 

mm 

32 

1.02 

1.02 

1.02 

mSm 

1.02 

1.02 

1.03 

1.02 

■BE 

40 

1.02 

1.02 

1.02 

1.03 

1.02 

1.02 

1.03 

1.03 

1.03 

48 

1.02 

1.02 

1.02 

1.03 

1.02 

1.02 

1.03 

1.03 

1.03 

56 

1.02 

1.02 

1.02 

1.03 

1.02 

1.02 

1.03 

1.03 

1.03 

64 

1.02 

1.02 

1.02 

1.03 

1.02 

1.02 

1.03 

1.03 

1.03 


Table 7. Orthotropy Parameters for [(0/90/±45) m ] s and [(0/90/+45 2 ) m ] s Laminates 


Number of 

Orthotropy parameter, 0 

laminate 

IM7 / 

Boron- 

S -glass- 

Kevlar 

AS4/ 

AS4/ 

Boron- 

IM7/ 


plies, 8m 

5260 

A1 

epoxy 

49-epoxy 

3502 

3501-6 

epoxy 

PETI-5 


8 



mm 


1.13 

1.13 

1.14 

1.14 

1.15 

16 



mSm 


1.07 

1.07 




24 

1.06 


1.04 


1.06 

1.06 

1.06 

1.06 


32 

1.05 


1.04 

1.06 

1.05 

1.05 

1.05 

1.05 


40 

1.04 


1.03 

1.05 

1.04 

1.04 

1.05 

1.05 


48 

1.04 


1.03 

1.05 

1.04 

1.04 

1.05 

1.04 


56 

1.04 


1.03 

1.04 

1.04 

1.04 

1.04 

1.04 


64 

1.04 


1.03 

1.04 

1.04 

1.03 

1.04 

1.04 



24 





































Table 8. Anisotropy Parameters for [(+<t>)2 m ]s Laminates 



Anisotropy par ameter, 

IM7/ 

5260 

Boron- 

A1 

S -glass- 
epoxy 

Kevlar 
49 -epoxy 

AS4/ 

3502 

AS4/ 

3501-6 

Boron- 

epoxy 

IM7/ 

PETI-5 


0 

LOO 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


5 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

mm, 

10 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

mm 

15 

1.00 

1.00 

1.00 

1.00 

1.00 

1.01 

1.01 

1.00 

m 

20 

1.00 

1.00 

1.00 

1.00 

1.00 

1.02 

1.01 

1.00 

mm 

25 

1.01 

1.00 

1.00 

1.00 

1.01 

1.03 

1.00 

1.00 

ail 

30 

1.02 

1.00 

1.00 

1.00 

1.01 

1.04 

1.00 

1.01 

all 

35 

1.03 

1.00 

1.01 

1.01 

1.02 

1.06 

1.00 

1.03 

apt 

40 

1.06 

1.00 

mm 

1.03 

1.04 

1.08 

1.01 

1.05 

mm 

45 

1.09 

1.00 

mm 

1.06 

1.07 

1.11 

1.05 

1.09 

aBB 

50 

1.13 

1.01 

1.06 

1.12 

1.11 

1.14 

1.11 

1.13 

mm 

55 

1.18 

1.01 

1.08 

1.18 

1.16 

1.17 

1.19 

1.19 

19 

60 

1.22 

mm 

1.09 

1.25 

1.20 

1.20 

1.27 

1.25 

1.33 

65 

1.25 

mm 

1.09 

1.29 

1.22 

1.21 

1.32 

1.29 

1.39 

70 

1.24 

1.02 

1.08 

1.29 

1.21 

1.19 

1.32 

1.29 

1.42 

75 

1.20 

1.01 

1.06 

1.24 

1.16 

1.15 

1.27 

1.24 

1.39 

80 

1.12 

1.01 

1.03 

1.15 

1.10 

1.09 

1.17 

1.15 

1.28 

85 

1.04 

1.00 

1.01 

1.05 

1.03 

1.03 

1.05 

1.05 

1.10 

90 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


25 


















Table 9. Anisotropy Parameters for | ±0 1-[- Laminates 


Fiber 

Anisotropy par ameter, 


IM7/ 

Boron- 

S -glass- 

Kevlar 

AS4/ 

AS4/ 

Boron- 

IM7/ 



5260 

A1 

epoxy 

49 -epoxy 

3502 

3501-6 

epoxy 

PETI-5 


0 


1.00 








5 


1.00 








10 

0.92 

1.00 

0.98 


0.94 

0.94 

0.90 

0.90 

0.85 

15 

0.88 

0.99 

0.96 

0.86 

0.90 

0.90 

0.86 

0.86 

0.80 

20 

0.85 

0.99 

0.94 

0.83 

0.87 

0.87 

0.83 

0.83 

0.77 

25 

0.83 

0.99 

0.93 

0.82 

0.85 

0.85 

0.82 

0.81 

0.76 

30 

0.82 

0.99 

0.92 

0.81 

0.84 

0.84 

0.82 

0.81 

0.75 

35 

0.82 

0.99 

0.92 

0.82 

0.84 

0.84 

0.82 

0.81 

0.75 

40 

0.83 

0.99 

0.93 

0.82 

0.85 

0.84 

0.83 

0.81 

0.75 

45 

0.84 

0.99 

0.94 

0.84 

0.86 

0.85 

0.85 

0.82 

0.76 

50 

0.86 

1.00 

0.95 

0.86 

0.88 

0.86 

0.87 

0.84 

0.77 

55 

0.88 

1.00 

0.96 

0.89 

0.90 

0.88 

0.90 

0.87 

0.78 

60 

0.91 

1.00 

0.98 

0.92 

0.93 

0.91 

0.94 

0.90 

0.80 

65 

0.94 

1.00 

0.99 

0.95 

0.96 

0.93 

0.97 

0.93 

0.84 

70 

0.97 

1.00 

0.99 

0.97 

0.98 

0.96 

0.99 

0.97 

0.89 

75 

0.99 

1.00 

1.00 

0.99 

0.99 

0.98 

1.00 

0.99 

0.94 

80 

1.00 

1.00 

1.00 

1.00 

1.00 

0.99 

1.00 

1.00 

0.98 

85 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

90 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


26 





















Table 10. Orthotropy Parameters for [0p/90g] T Laminates 


P x 100% 
p + q 

Orthotropy parameter, 0 

IM7/ 

5260 

Boron- 

A1 

S -glass- 
epoxy 

Kevlar 

49-epoxy 

AS4/ 

3502 

AS4/ 

3501-6 

Boron- 

epoxy 

IM7/ 

PETI-5 


0 


MEM 

0.69 

0.52 

0.55 

■j&B 

0.55 

0.49 

0.34 

5 



0.77 

0.68 

0.69 


0.69 

0.66 

0.63 

10 

mm 

m m 

0.82 

0.76 

0.77 

mm 

0.77 

0.75 

0.73 

15 

0.81 

0.95 

0.86 

0.81 

0.82 

0.81 

0.82 

0.81 

0.79 

20 

0.85 

0.97 

0.89 

0.86 

0.86 

0.85 

0.86 

0.85 

0.84 

25 

0.89 

0.98 

0.92 

0.89 

0.89 

0.89 

0.89 

0.89 

0.88 

30 

0.92 

0.98 

0.94 

0.92 

0.92 

0.91 

0.92 

0.91 

0.91 

35 

0.94 

0.99 

0.96 

0.94 

0.94 

0.94 

0.94 

0.94 

0.93 

40 

0.96 

1.00 

0.97 

0.96 

0.96 

0.96 

0.96 

0.96 

0.96 

45 

0.98 

1.00 

0.99 

0.98 

0.98 

0.98 

0.98 

0.98 

0.98 

50 

1.00 

1.01 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

55 

1.02 

1.01 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

1.03 

60 

1.05 

1.02 

1.04 

1.05 

1.05 

1.05 

1.05 

1.05 

1.06 

65 

1.09 

1.03 

1.06 

1.08 

1.08 

1.09 

1.08 

1.09 

1.10 

70 

1.13 

1.04 

1.09 

1.13 

1.12 

1.13 

1.12 

1.13 

1.15 

75 

1.19 

1.05 

1.12 

1.18 

1.17 

1.19 

1.17 

1.19 

1.22 

80 

1.26 

1.06 

1.17 

1.25 

1.24 

1.26 

1.24 

1.27 

1.31 

85 

1.35 

1.07 

1.22 

1.35 

1.33 

1.36 

1.33 

1.37 

1.43 

90 

1.48 

1.09 

1.28 

1.47 

1.44 

1.49 

1.44 

1.50 

1.61 

95 

1.67 

1.11 

1.36 

1.65 

1.60 

1.67 

1.59 

1.70 

1.93 

100 

1.98 

1.13 

1.45 

1.93 

1.84 

1.98 

1.82 

2.05 

2.93 
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Table 1 1 . Anisotropy Parameters for [0p/90g] T Laminates 


P x 100% 
p + q 

Anisotropy par ameter, 

IM7/ 

5260 

Boron- 

A1 

S -glass- 
epoxy 

Kevlar 49- 
epoxy 

AS4/ 

3502 

AS4/ 

3501-6 

Boron- 

epoxy 

IM7/ 

PETI-5 


0 


1.00 

1.00 







5 


1.00 

0.98 







10 


0.99 

0.95 


0.87 




0.64 

15 


0.99 

0.93 


0.83 


0.84 

0.78 

0.59 

20 


0.98 

0.91 

0.78 

0.81 


0.81 

0.75 

0.57 

25 

0.75 

0.98 

0.90 

0.76 

0.79 

0.75 

0.80 

0.73 

0.57 

30 

0.74 

0.98 

0.89 

0.75 

0.78 

0.74 

0.79 

0.73 

0.59 

35 

0.75 

0.98 

0.88 

0.75 

0.78 

0.74 

0.78 

0.73 

0.61 

40 

0.75 

0.97 

0.88 

0.76 

0.78 

0.75 

0.79 

0.74 

0.65 

45 

0.77 

0.97 

0.88 

0.77 

0.79 

0.77 

0.80 

0.76 

0.69 

50 

0.79 

0.97 

0.89 

0.79 

0.81 

0.79 

0.82 

0.78 

0.73 

55 

0.82 

0.97 

0.90 

0.82 

0.83 

0.82 

0.84 

0.81 

0.77 

60 

0.85 

0.98 

0.91 

0.85 

0.86 

0.85 

0.86 

0.84 

0.82 

65 

0.88 

0.98 

0.92 

0.88 

0.89 

0.88 

0.89 

0.88 

0.86 

70 

0.91 

0.98 

0.94 

0.91 

0.91 

0.91 

0.92 

0.91 

0.90 

75 

0.94 

0.98 

0.95 

0.93 

0.94 

0.94 

0.94 

0.94 

0.93 

80 

0.96 

0.99 

0.97 

0.96 

0.96 

0.96 

0.96 

0.96 

0.96 

85 

0.98 

0.99 

0.98 

0.97 

0.98 

0.98 

0.98 

0.98 

0.98 

90 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

95 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

100 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 
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Table 12. Orthotropy Parameters for 1 7() |v /() l| |-] Laminates 


P x 100% 
p + q 

Orthotropy parameter, 0 

IM7/ 

5260 

Boron- 

A1 

S -glass- 
epoxy 

Kevlar 49- 
epoxy 

AS4/ 

3502 

AS4/ 

3501-6 

Boron- 

epoxy 

IM7/ 

PETI-5 


0 

MEM 

1.13 

WSM 

1.93 

mm 

mm 

mm 

2.05 

2.93 

5 

IBS 

1.11 

Bn 

1.69 

mm 

m 

IP 

1.75 

2.02 

10 

mm 

1.09 

US 

1.52 

wm 

mm 

mm 

1.56 

1.70 

15 

1.41 

1.08 

1.24 

1.41 

1.38 

1.41 

1.38 

1.43 

1.52 

20 

1.32 

1.07 

1.20 

1.31 

1.30 

1.32 

1.30 

1.33 

1.39 

25 

1.25 

1.05 

1.16 

1.24 

1.23 

1.25 

1.23 

1.25 

1.29 

30 

1.19 

1.05 

1.13 

1.19 

1.18 

1.19 

1.18 

1.20 

1.22 

35 

1.15 

1.04 

1.10 

1.15 

1.14 

1.15 

1.14 

1.15 

1.17 

40 

1.11 

1.03 

1.08 

1.11 

1.11 

1.11 

1.11 

1.12 

1.13 

45 

1.09 

1.03 

1.06 

1.09 

1.08 

1.09 

1.08 

1.09 

1.10 

50 

1.06 

1.03 

1.05 

1.06 

1.06 

1.06 

1.06 

1.06 

1.07 

55 

1.04 

1.02 

1.03 

1.04 

1.04 

1.04 

1.04 

1.04 

1.05 

60 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

65 

1.00 

1.01 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

70 

0.98 

1.01 

0.99 

0.98 

0.98 

0.98 

0.98 

0.98 

0.97 

75 

0.95 

1.00 

0.97 

0.95 

0.95 

0.95 

0.96 

0.95 

0.95 

80 

0.92 

0.99 

0.94 

0.92 

0.92 

0.92 

0.92 

0.92 

0.91 

85 

0.88 

0.98 

0.91 

0.88 

0.88 

0.88 

0.89 

0.88 

0.86 

90 

0.82 

0.97 

0.88 

0.83 

0.83 

0.82 

0.83 

0.82 

0.80 

95 

0.74 

0.95 

0.83 

0.75 

0.76 

0.74 

0.76 

0.74 

0.71 

100 

0.60 

0.93 

0.75 

0.61 

0.63 

0.60 

0.62 

0.59 

0.51 
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Table 13. Anisotropy Parameters for [ 7() |v /() l| |-] Laminates 


P x 100% 
p + q 

Anisotropy par ameter, 

IM7/ 

5260 

Boron- 

A1 

S -glass- 
epoxy 

Kevlar 49- 
epoxy 

AS4/ 

3502 

AS4/ 

3501-6 

Boron- 

epoxy 

IM7/ 

PETI-5 


0 

1.00 


1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


5 

1.02 

1.01 

1.01 

1.02 

1.02 

1.02 

1.02 

1.03 

wEm 

10 

1.04 

1.01 

mgm 

1.05 

1.04 

1.04 

1.05 

1.05 

mm 

15 

1.06 

1.02 

m S 

1.07 

1.05 

1.05 

1.07 

1.08 

Up 

20 

1.07 

1.02 

1.03 

1.09 

1.06 

1.05 

1.09 

1.09 

m 

25 

1.07 

1.02 

1.03 

1.10 

1.06 

1.05 

1.10 

1.10 

■P 

30 

1.07 

1.03 

1.03 

1.10 

1.06 

1.04 

1.11 

1.09 

mm 

35 

1.06 

1.03 

1.03 

1.10 

1.05 

1.03 

1.11 

1.09 

1.18 

40 

1.05 

1.03 

1.03 

1.09 

1.04 

1.01 

1.12 

1.08 

1.16 

45 

1.03 

1.04 

1.03 

1.09 

1.04 

1.00 

1.12 

1.06 

1.12 

50 

1.02 

1.04 

1.03 

1.08 

1.03 

0.98 

1.12 

1.05 

1.08 

55 

1.01 

1.04 

1.03 

1.08 

1.02 

0.97 

1.13 

1.04 

1.05 

60 

1.01 

1.04 

1.03 

1.08 

1.02 

0.96 

1.14 

1.04 

1.01 

65 

1.01 

1.04 

1.04 

1.09 

1.03 

0.96 

1.15 

1.04 

0.98 

70 

1.02 

1.04 

1.05 

1.11 

1.04 

0.97 

1.18 

1.05 

0.96 

75 

1.04 

1.04 

1.06 

1.13 

1.06 

0.98 

1.20 

1.07 

0.96 

80 

1.07 

1.04 

1.07 

1.16 

1.09 

1.00 

1.24 

1.10 

0.97 

85 

1.11 

1.04 

1.08 

1.20 

1.12 

1.04 

1.28 

1.15 

1.01 

90 

1.16 

1.03 

1.09 

1.25 

1.16 

1.09 

1.32 

1.20 

1.08 

95 

1.22 

1.03 

1.10 

1.31 

1.20 

1.15 

1.36 

1.27 

1.21 

100 

1.24 

1.02 

1.08 

1.29 

1.21 

1.19 

1.32 

1.29 

1.42 


30 




















Nondimensional 


90%-decay length, 


d 

^Rh 


8 = 0.1 



Orthotropy parameter, 0 = 


12 AA, 
(A„A 22 -A’ 2 )h s 


1/4 


Figure 1 . Nondimensional 90-percent-decay length for symmetrically and unsymmetrically laminated cylinders, 
as a function of laminate orthotropy. 



Anisotropy parameter, 


Figure 2. Nondimensional 90-percent-decay length for symmetrically and unsymmetrically laminated cylinders, 
as a function of laminate anisotropy. 
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Figure 3. Effect of lamina material properties on nondimensional orthotropy parameter for single-ply, homoge- 
neous, specially orthotropic laminates (0.2 < v 12 < 0.35). 



Fiber angle deg 


Figure 4. Effect of lamina material properties on nondimensional orthotropy parameter for [(±(|)) in ] s , t(+(j>)2 m ls’ an d 
[(±(|)) m ] T laminates (m = 1, 2, ...). 
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0 8 16 24 32 40 48 

Number of laminate plies, 8m 

Figure 5. Nondimensional orthotropy parameter for typical laminates made of IM7/5260 material. 



Number of laminate plies, 8m 

Figure 6. Effect of lamina material properties on nondimensional orthotropy parameter for [(0 2 /±45) m ] s laminates. 
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P-1 00/3502 IM7/PETI-5, Kevlar 49-epoxy, 
and f 


Orthotropy 

parameter, 


0 = 


Boron-epc 

AS4/3502 


l ! ! 0.4 


tS4/3501 -6 
j IM7/5260 


Kevlar 49-epoxy 


12A.D,, 

1/4 

1.0 

I 1 , | 1 


— 


(A„A 22 - A, 2 )h 2 


"|| ff% 

i S-giass-epoxy 

It 

ilass-c 


Boron- 


P-1 00/3502 


[(±45/0/90) m ] s and [(+45 2 /0/90) m ] s , 
- [(0/90/±45) m ] s and [(0/90/+45 2 ) m ] s 


•8| i L 


16 24 32 

Number of laminate plies, 8m 


40 


48 


Figure 7. Effect of lamina material properties on nondimensional orthotropy parameter for quasi-isotropic 
laminates and similar unbalanced laminates. 



Figure 8. Effect of lamina material properties on nondimensional anisotropy parameter for [(+<j>) 2 m ] s laminates 
(m = 1, 2, ...). 
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Figure 9. Effect of lamina material properties on nondime nsional first-order correction factor for [(+<t>)2m]s 
laminates (m =1,2, ...). 


Nondimensional 


90%-decay length, 


d 

*/Rh 


£= 0.1 



Figure 10. Nondimensional 90-percent-decay length for [(+0)2 m l s laminates made of IM7/5260 graphite - 

1 


bismalehnide and P-100/3502 pitch-epoxy material I — = m = 1, 2, 

VR 20 


35 




Anisotropy 

parameter, 



Number of laminate plies, 8m 


Figure 1 1 . Effect of lamina material properties on nondimensional anisotropy parameter for [ (0/90/+45 2 )i n ] s and 
[(+45 2 /0/90) m ] s laminates. 


Anisotropy 

parameter, 



Fiber angle 4>, deg 


Figure 12. Effect of lamina material properties on nondimensional anisotropy parameter for | ■ 0]j laminates. 
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1.2 


1.0 


.8 


Anisotropy 
parameter, .6 


.4 


.2 


0 

0 15 30 45 60 75 90 

Fiber angle <j>, deg 

Figure 13. Effect of stacking sequence number on nondimensional anisotropy parameter for [(±<])) m ] T laminates 
made of P-100/3502 pitch-epoxy material. 



Fiber angle 4>, deg 

Figure 14. Effect of lamina material properties on nondimensional first-order correction factor for | laminates. 
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Figure 15. Effect of stacking sequence number on nondimensional first -order correction factor for [(±(f>) m ] T 
laminates made of P-100/3502 pitch-epoxy material. 



Figure 16. Nondimensional 90-percent-decay length for [d=c]>l T laminates made of IM7/5260 graphite -bismaleimide 
and P-100/3502 pitch-epoxy material 
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Percentage of 0° plies, p + q x 100% 


Figure 17. Effect of lamina material properties on nondimensional orthotropy parameter for [0p/90„] T laminates. 


Anisotropy 

parameter, 

do 



Percentage of 0° plies, p^q x 100% 


Figure 18. Effect of lamina material properties on nondimensional anisotropy parameter for [0p/90 q ] T laminates. 
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Nondimensional 
30%-decay len 
d 

^Rh „ = „, 



Figure 19. Nondimensional 90-percent-decay length for [0p/90 q ] T laminates made of IM7/5260 graphite - 

bismaleimide and P-100/3502 pitch-epoxy material [ ^ 

V R 20 J 



Figure 20. Effect of lamina material properties on nondimensional orthotropy parameter for [70p/0 q ] T laminates. 
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Figure 21. Effect of lamina material properties on nondimensional anisotropy parameter for [70 |) /0 (| ]'p laminates. 


Nondimensional 
90%-decay length, 

d 

^Rh e=01 



Figure 22. Nondimensional 90-percent-decay length for [70^0 ] T 

(h = 1_) 

VR 20/ 


bismalehnide and P-100/3502 pitch-epoxy material 


laminates made of IM7/5260 graphite- 
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Appendix A 


Sanders-Koiter Equations 


The linear Sanders-Koiter shell equations (refs. 16 and 17) are presented in this appendix for a 
right-circular cylinder with a radius that is given by R. For these equations, x and 0 denote the axial 
and circumferential coordinates, respectively. First, the equilibrium equations are presented; then the 
kinematic equations and the constitutive equations are presented. Last, the boundary conditions are 
given for a complete right-circular cylinder at an edge that is given by a constant value of the axial 
coordinate x. 


Equilibrium Equations 

The equilibrium equations are given in a form similar to those found in reference 22; that is, 


3N X i 3N X0 c 2 3M X0 
3x + R 30 2R 2 30 +qx_ 


(Al) 


3N X0 i 3N 0 Ci c 2 3M X0 

“I _ 1 O ( ) H + Q Q — 0 

3x R 30 R 0 2R 3x 0 


(A2) 


0 Q X | l 0 Qe 

~dx + R 30 


N e 

R 


+ q n = 0 


(A3) 


3M X i 3M X0 
3x + R 30 

3M X0 i 3M 0 
3x + R 30 


-Q x =o 

-Qe=0 


(A4) 

(A5) 


where N x , N 0 , and N x0 are the membrane stress resultants; Q x and Q 0 are the transverse shear- 
stress resultants; M x , M 0 , and M x0 are the bending stress resultants; q x , q 0 , and q n are the applied 
surface tractions; and c 1 and c 2 are constants that identify the equations of other shell theories that 
are considered herein. In particular, the Sanders-Koiter equations are given by c^c^l, and the 
Love-Kirchhoff equations are given by c ^ = 1 and c 2 = 0 . Donnell’ s equations are given by c ! = c 2 = 0 . 
This convention is used throughout the present study. 

Kinematic Equations 

The kinematic equations are given by 


P o 3u 
x 3x 


(A6) 
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0 _ 1 3v w 
80 -R90 + R 


n _3v 1 3ll 

< ™”dx + R 30 


3 o dw 


Px=- 


dx 


nO C 1 1 dw 

p « = r v -r59 


ft=- 


c 2[dv 1 3u) 


2 3 x Rae 




o 

.o ” r 'X 


K x = _ 
x dx 


a 2 w 


o _ l 3 Pe cj 3 V ^ a 2 w 
K0_ r ae R - 


» 2 ae r 2 ae 2 


i / 

O _ 1 Kx , oo 

Kx0 'r ^ +Pn 


J 


aPe. 

3x 


2 a 2 w 


1 


r axae r 


Cl + 


)3v 

,3x 


z 2 3u 


2R 


2 30 


(A7) 

(A8) 

(A9) 

(A10) 

(All) 

(A12) 

(A13) 

(A14) 


where u, v, and w are the axial, circumferential, and radial displacements of a point of the shell middle 
surface; 8°, e 0 , and y"q are the membrane strains; |3°, (3 0 , and (3° are the rotations; and K x , K 0 , and 
K° 0 are the bending strains. The displacement w is positive when it is outward from the cylinder 
reference surface. 


Constitutive Equations 

The isothermal constitutive equations are given in matrix form by 


N A 


A 11 

A 12 

A 16 

Bn 

B 12 

B 16 

N e 


A 12 

A 22 

A 26 

B 12 

B 22 

B 26 

N X0 

\- 

A 16 

A 26 

A 66 

B 16 

B 26 

b 66 

M x 

[- 

Bn 

B 12 

B 16 

Dll 

D 12 

Die 

M e 


B 12 

B 22 

B 26 

D 12 

D 22 

D 26 

V M X0 j 


B 16 

B 26 

B 6 6 

Die 

D 26 

d 66 


( ~o\ 


£e 

7x0 




4 

j 


(A15) 
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where the subscripted A, B, and D terms of the matrix are the stiffnesses of laminated-composite 
shells that are obtained from the Love-Kirchhoff shell theory. Moreover, the constitutive terms in equa- 
tion (A15) are identical to those for laminated-composite plates that are given in reference 18, page 198. 

Boundary Conditions 

The boundary conditions for an edge that is defined by a constant value of the axial coordinate x 
are given by 


N x - N x (0) or u = u(e) 


(A 16) 


N. 


x0 ' 


R 


Ci + 


M x0 =T(0) or v = v(0) 


(A- 17) 


1 dM x 0 _ 

Q x +— -^-=V(0) or w = w(0) (A 18) 

M x = M x (0) or P° = p(0) (A19) 

where u(0) , v(0) , and w(0) are applied edge displacements; P(0) is an applied edge rotation; and N x (0) , 
T(0), V(0) , and M x (0) are applied edge loads. 
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Appendix B 

Equations for Axisymmetry 


The linear Sanders-Koiter shell equations that are presented in appendix A for a right-circular cylin- 
der with a radius R are specialized in this appendix for the case of axisymmetric behavior. For these 
equations, x and 0 denote the axial and circumferential coordinates, respectively. The specialization to 
axial symmetry is conducted by eliminating all terms in the equations of appendix A that are differenti- 
ated with respect to the circumferential coordinate 0. First, the equilibrium equations, the kinematic 
equations, and the constitutive equations are presented. Then, the boundary conditions are given for a 
complete right-circular cylinder at an edge that is given by a constant value of the axial coordinate x. 
Last, the axisymmetric equations are manipulated into a single ordinary differential equation in terms of 
the radial displacement w(x). 

Equilibrium Equations 

The equilibrium equations for axisymmetric behavior are given by 

dN x 

-+q x (x) = 0 (Bl) 


dN x0 C 1 _ c 2 dM x0 

H Qn + 

dx R B 2R dx 


+ q e (x) = 0 


(B2) 


dQx 

dx 


N e 

R 


+ q n (x)=0 


(B3) 


dMx 

dx 


Q x =o 


(B4) 


dM 


X0 


dx 


-Q e = 0 


(B5) 


where the membrane stress resultants N x , N 0 , and N x0 ; the transverse shear- stress resultants Q x and 
Q 0 ; the bending stress resultants M x , M 0 , and M x0 ; and the applied surface tractions q x , q 0 , and q n 
are functions of only the axial coordinate x. 

Kinematic Equations 

The kinematic equations are given by 


o du 

x dx 


(B6) 
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80 = — 
9 R 


(B7) 


y£ 0 =— 
fx0 dx 


(B8) 


PS=-^ 

Hx dx 


(B9) 


C 1 

(3 0 = — v 
R 


(BIO) 


50 c 2 dv 


Pn = 


2 dx 


(B 11) 


d P* 

K° = 

x dx 


d 2 w 

dx 2 


(B12) 


«e=o 


(B13) 


1 d(3 0 i 

k° 0 =-|3° + — =- 
x0 R Hn dx R 


1 \dv 
Ci H — Co 
1 2 2 dx 


(B14) 


where the middle-surface displacements u, v, and w; the membrane strains 8°, 8 0 , and y£ 0 ; the 
rotations (3°, (3 0 , and (3°; and the bending strains K°, K 0 , and K° 0 are functions of only the axial 
coordinate x. 


Constitutive Equations 

The isothermal constitutive equations reduce to 


N A 


A 11 

A 12 

A 16 

Bn 

B 12 

B 16 

N e 


A 12 

A 22 

A 26 

B 12 

B 22 

B 26 

N X0 

\- 

A 16 

A 26 

A 66 

B 16 

B 26 

b 66 

M x 

f- 

B it 

B 12 

B 16 

Dll 

D 12 

Die 

M e 


B 12 

B 22 

B 26 

D 12 

D 22 

D 26 

J 


B 16 

B 26 

B 6 6 

Die 

D 26 

d 66 


f e o\ 

C X 

8 e 

Yxe 


j 


(B15) 
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where the subscripted A, B, and D terms of the matrix are the usual constitutive terms of classical 
Love-Kirchhoff-type laminated composite shell theory or classical laminated plate theory (e.g., see 
p. 198 of ref. 18). 

Boundary Conditions 

The boundary conditions for an edge that is defined by a constant value of the axial coordinate x 
are given by 


N X =N X or u = u 


(B16) 


N X0 + 


R 


c l+- c 2 M xe= T or v = v 


(B17) 


Q = V or w = w 


(B 18 ) 


M X = M X or Px = P (B19) 

where the applied edge displacements u, v, and w; the applied edge rotation (3; and the applied edge 
loads N x , T, V, and M x are all constants. 


Bending Boundary-Layer Equation 

The bending boundary-layer equation is obtained by first noting that integration of equation (Bl) 
yields 


N 


X 



q x dx + C=N(x) 


(B20) 


where C is a constant of integration that is determined from the boundary condition given by equa- 
tion (B16). Next, equations (B2) and (B5) are combined to get 


dN x9 

dx 


R 


1 

H — 

2 


c 2 I 


dM 


X0 


dx 


+ q e (x)=0 


(B21) 


For convenience, the parameter 


1 

M’ =c i + 2 c 


2 


(B22) 


is introduced such that the Sanders-Koiter equations are given by |l = — and the Love-Kirchhoff equa- 
tions are given by |J, = 1 . Donnell’s equations are given by ll = 0 . Similarly, the function 
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(B23) 


T(x) = N xe+ |M xe 

is introduced so that equation (B21) becomes 

dT . , . 

^ + q e » = ° 

and the corresponding boundary condition given by equation (B17) becomes 

T = T or v =v 


Integration of equation (B24) yields 

T=-Jq 0 dx + C= T(x) 


(B24) 


(B25) 


(B26) 


where C is a constant of integration that is determined from the boundary condition given by equa- 
tion (B25). Next, equations (B3) and (B4) are combined to give 


d 2 M„ 


dx- 


N e 

R 


+ q n (x) = 0 


(B27) 


The next step in the analysis is the simplification of the constitutive equations. First, by using equa- 
tions (B 8 ) and (B22), equation (B14) is expressed as 


By using equations (B23) and (B28), the constitutive equations are expressed as 


N I 


A 11 

A 12 

A 16 

B 11 


fe°) 

C X 

N 0 

V- 

A 12 

A 22 

A 26 

B 12 

j 


f 1 


A 16 

A 26 

A 66 

B 16 

| 

^0 

w 


B 11 

B 12 

B 16 

Dll 


K° 

\ X J 


M 0 - B 12 e x + B 22 e e + B 26^0 + D 12 k x 


(B28) 


(B29) 


(B30) 
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where 


A 16 = A 16 "*"M' 


f h} 

1 R J 


B 


16 


(B31) 


A 26 = A 26 ■*" M' 


fu\ 


B 


26 




(B32) 


A 66 - A 66 + 


fh\ 


B 


v R y 


66 2i 

— +|1 1 




2 D, 


l R ; 


66 

2~ 


(B33) 


B i6- B i6 + ^ 


r h 

i R ; 


D 


16 


B 26 — B 26 + 


fh\ 


v R j 


D 


26 


(B34) 


(B35) 


The motivation for writing the constitutive equations in this form is that the matrix equation given by 
equation (B29) is the only part of the full constitutive equations that appears in the strain-energy density 
function, which is used in the present paper to determine the corresponding positive-definiteness condi- 
tions. With these simplified constitutive equations and equations (B6)-(B8) and (B12), equation (B20) 
is expressed as 


du 


w 


dv 


d 2 w 


Aii di +Ai2 R +Ai6 *r Bu d ^‘“ N<x> ~° 


(B36) 


and equation (B26) is expressed as 


— du t w t dv - d w - 
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dx 


+ A 26 ^ +A 66 “^“B 


dx 


16 


dx 


-T(x)=0 


(B37) 


du dv 

Equations (B36) and (B37) are then solved for — and — to get 


du 

dx 


i t w . i 

A 66 N ( x ) ““ A 16 T ( x ) + ^ A 16 A 26 “ A 12 A 66 Jr + ( A 66 B 11 “ A 16 B 16j 


td 2 w 


dx z 


A 11 A 66“ A 16 


(B38) 


dv 

dx 


9 

- — — / — — \ w / — — \d^w 

A llT(x)- A 16 N(x) + |A 12 A 16 - Aj^gJ— + |A n B 16 - A 16 B n j - 


dx 


— — 2 

A 11 A 66“ A 16 
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(B39) 



Equation (B39) indicates that the circumferential displacement v(x) becomes uncoupled from the axial 
displacement u(x) and the radial displacement w(x) when A 16 = A 2 6 = B 16 = 0, which implies that 
Ai6 = A 26 = Bi6 = Di6 = 0. In addition, the constitutive equation (B29) indicates that N x , N 0 , and 
M x become uncoupled from the torsional, shear strain 7x0 when A ] ( , = A 2f) = B i 6 =0 and that T, 
which is defined by equation (B23), becomes uncoupled from 8°, e 0 , and K x . Furthemiore, equa- 
tion (B 30) indicates that M 0 becomes uncoupled from 7x0 when B 26 = 0, which implies 

B 26 = D 26 = 0 - 

Next, equations (B38) and (B39) are then substituted into equations (B6) and (B8), and the resulting 
expressions for 8° and y^ 0 , along with equations (B7) and (B12), are substituted into the constitutive 
equation (B29). This action converts the strains and stress results in equation (B29) into functions of the 
radial displacement w(x). Substituting the expressions for N 0 and M x into equation (B27) yields the 
bending boundary-layer equation that is given by 

^ d 4 w _ d 2 w _ _ , . 

Cq — T + C 2 — - + C 3 w = C 4 (x) (B40) 

dx 4 dx 2 


The constant coefficients are given by 


Ci -D u 


1 - 


A 66 B 11 + A 11 B 16 _2A 16 B 11 B 16 


A 11 A 66~ A 16) D U 


(B41) 


Co =— : 


R 


B 12 + 


A 16 A 26“ A 12 A 66 B 11 + A 12 A 16~ A 11 A 26 B 16 


A 1 1 A 66 A 1 6 


(B42) 


c 3 = 


A 11 A 22 — A 12 A 66 “ A 11 A 26 _ A 22 A 16 +2A 12 A 16 A 26 


R- 


A 1 1 A 66 A 16 


(B43) 


The function C 4 (x) is given by 


C 4 (x) = q n (x) + 


A 16 A 26 “ A 12 A 66 N ( x ) + A 12 A 16 “ A 11 A 26 T ( x ) 




R l A 11 A 66“ A 16 


d 2 N 


B 1 1 A 66 A 1 6 B 16 2 + A 11 B 16 _B 11 A 16| 2 

dx 2 1 >dx 2 


) d 2 T 


— — 2 

A 1 1 A 66 A 16 


(B44) 
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These expressions are simplifed further by introducing the following expressions: 


a 12“ 


A 16 A 26 _A 12 A 66 


A 11 A 22 _ A 12 ) A 66 _ A 11 A 26 “ A 22 A 16 + 2A 12 A 16 A 26 


(B45) 


a 22 — 


A 11 A 66“ A 16 


A 1 1 A 22 ~ A 12 A 66 A 1 1 A 26 A 22 A 1 6 + 2A 1 2 A 1 6 A 26 


(B46) 


a 26 — 


A 1 2 A 16 A 11 A 26 

A 1 1 A 22 “ A 12 ) A 66 _ A 1 1 A 26 ~ A 22 A 16 + 2A 12 A 1 6 A 26 


(B47) 


t>21~- a 12 B ll +a 22 B 12 + a 26 B 


16 


(B48) 


0 = 1 - 


A 66 B 11 +A 11 B 16 _2A 16 B 11 B 16 


A 1 1 A 66 ' A 16| l ^11 




(B49) 


By using equations (B45)-(B49), equations (B41)-(B43) are expressed as 


Ci=D n e 


(B50) 


Co = 


2 h 

R a 


21 


22 


(B51) 


Co = 


R 2 a 


22 


(B52) 


Similarly, for the case where the second derivatives of N(x) and T(x) are zero valued, equation (B44) 
becomes 


C 4( x ) = c ln( x ) + 


a 12 N(x) + a 26 T(x) 
Ba 22 


(B53) 


The desired fomr of the bending boundary-layer equation is obtained by dividing equation (B40) by C p 
that is, 
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(B54) 


4 2 

d w d“w 

— 7 +4S — T 

dx 4 dx" 


+ 4Qw = P(x) 


where the constants S and Q are given by 



4C 1 2Ra 22 D n e 


(B55) 


Q = 


C3 

4Ci 


1 

4R 2 a 22 D 11 e 


The function P(x) is given by 


C 4 (x) q n (x) a 12 N(x) + a 26 T(x) 

P(x) = = + — 

C-i D n e Ra 22 Djje 


(B56) 


(B57) 


for the special case when the second derivatives of N(x) and T(x) are zero valued. The quantity D n e 
that appears in equations (B54)-(B56) is sometimes referred to, in some contexts, as a reduced bending 
stiffness (ref. 12). 
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Appendix C 

Anisotropy-Factor Equations 

The first-order approximation of the anisotropy factor /t that is used herein is given by 


Q 





(Cl) 


where is the value of equation (20b), with |A = 0 in the terms with the overbars. This expression is 
given by 




A 11 A 22“ A 12 


L 11 


a 22 e 0 


Y l4 f 

) 


- 1/2 


1- 


} 21 


V 




a 22 D ll e 0 


(C2) 


which is the anisotropy factor that corresponds to Donnell’s equations. The terms a 22 , b 2 p and e () are 
given by 


a 12“ 


A 16 A 26“ A 12 A 66 


A 11 A 22“~ A 12 ) A 66 — A 11 A 26~ A 22 A 16 + 2A 12 A 16 A 26 


(C3) 


a 22 “ 7 


A 1 1 A 66 A 1 6 


A 1 1 A 22 A 12 A 66“ A 11 A 26 _ A 22 A 16 +2A 12 A 16 A 26 


a 26 — 


A 12 A 16 A 1 1 A 26 

A 11 A 22 — A 12 A 66“ A 11 A 26 _ A 22 A 16 +2A 12 A 16 A 26 


b 21 


a 12 B ll +a 22 B 12 + a 26 B 16 


e Q - 1 - 


A 66 B 11 + A 11 B 16 — 2A 16 B 11 B 16 
A 11 A 66 -A 16 D 11 


(C4) 


(C5) 


(C6) 


(C7) 
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The term is a first-order correction to the results that correspond to Donnell’s equations and is 
given by 


<?i = 


\j a 22 D 1 l e 0 ( a 22 e 1 +a 122 e o] +2a 22^ b 121 e 0“ b 21 e l j “ 2a 122 b 21 e 0 

4a 22 e 0 (y^ a 22 D ll e 0 — b 2lJ 


(C8) 


where 


2 I A 16 B 11“ A 11 B 16 


e, = 




D 


16 


1 


A 11 A 66~ A 16 


— B i6 B n A 66 + A n B 66 + A 16 A n B 66 +B 16 


A 


( 


A 11 A 66~ A 16 D ll b 


(C9) 


and 


a 1 12 _ 


B 1 6 f 1 6 + B 26*26 + B 66*66 


A 11 A 22~ A 12 A 66^ A 11 A 26 A 22 A 16 +2A 12 A 16 A 26 


(CIO) 


a 122 2 


B 1 6§ 16 '* B 26g26 + B 66§66 


A 1 1 A 22 ^ A 12 A 66 ~ A 1 1 A 26 “ A 22 A 1 6 + 2A 1 2 A 1 6 A 26 


h 


(Cll) 


a 126 : 


B 16 b 16 + B 26 b 26 + B 66 b 66 


A 1 1 A 22 “ A 12 A 66 ~ A 1 1 A 26 “ A 22 A 1 6 + 2A 12 A 1 6 A 26 


h 


(Cl 2) 


b 121 - “ 


a 112Bll+ a l22 B l2 + a 126Bl6 + 


a 26 D 16 

h 


(Cl 3) 


f 16 - A 26| A 16 A 22~ A 11 A 26 


+ A, 


66 


A 26 A 11 A 22 + A 12 _ 2A 12 A 2 2A 
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(C14) 
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f 26 - A 16( A 26 A 11~ A 22 A 16, + A 66 A ie( A ll A 22 + A 12^ “ 2A 11 A 12 A 26 (C15) 

f 66 =2 ( A 12 A 16~ A 11 A 26)( A 16 A 22~ A 12 A 2e) (C16) 

§16 = ( A 12 A 16 - A 11 A 26^ A 16 A 26 — A 12 A 66^ (C17) 

§26 = | A 12 A 16 _ A ll A 26j[ A ll A 66“ A 16j (Cl 8) 

§66 = | A 12 A 16 — A 1 1 A 26^ (Cl 9) 


h 16 ~ A 12 A 66^ A 12 A 66 — A 12^ + A 1 1 A 26^ A 12 A 26 — A 16 A 22^ + A 16 A 22^ A 12 A 16 “ A ll A 2ej (C20) 
h 26 = ~ A 1 1 A 66^ A 1 1 A 22 — A 12 j + A 16 A 11 A 22 ^ 2A 12^ + A 1 1 A 26^ 2A 12 A 16 “ A 11 A 26^ (C21) 




A 11 A 22 



(C22) 


Special Cases for /* 0 and 

Simplifications to and are presented below for unbalanced and balanced, symmetric lami- 
nates and for balanced, unsymmetric laminates that include the subclasses of general antisymmetric 
laminates, antisymmetric cross-ply laminates, and antisymmetric angle-ply laminates. 


Unbalanced and balanced symmetric laminates 

For unbalanced, symmetric laminates, A 16 ^0, A 2 6 : ? t 0, and B^ =B 1 2 = B 22 = B 16 = B 26 = 
B 66 = 0. For this special case, 
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which agrees with the corresponding equations given by Reuter (ref. 4), and 


e,= 


D 1 6 A 12 A 16 _ A 11 A 26 


2h 


A 11 A 22~ A 12 


A 66 _ A 1 1 A 26 “ A 22 A 16 + 2A 12 A 16 A 26 


1 / 2 , 


D lll A 11 A 66“ A 16 


1/2 


(C24) 


For balanced, symmetric laminates, A 16 = A 2 g = 0 in addition to the subscripted B-matrix constitutive 
terms. For this special case, /t ^ - 1 and ^ =0. 

Balanced , unsymmetric laminates 

For balanced, unsymmetric laminates, ^\6 = ^26 = 0, which yields the following simplified 
expressions: 


e \ 

° A 11 D 11 A 66 D 11 

A 11 B 12~ A 12 B 11 
\J A 1 1 D 1 l( A l 1 A 22 “ A 12) 

^o = \Ao( 1+a ) 1/2 


(C25) 


(C26) 


(C27) 


ej = -2 — 

1 A 66h 


f. 


D 16 B 11 B 16 B 16 B 66 


D 11 A 11 D 11 A 66 D 11 


V 


(C28) 


e,= 


i 


4(1 +A) 


(1 +2A) + 


2B 16 A 1 1 B 26 A 1 2 B 1 6 


° A 66 h \/ A 11 D 11 A 11 A 22 — A 12 ) e 0 


(C29) 


For the subclass of balanced, antisymmetric laminates, D 16 = D 2 g = 0, in addition to the shear- 
extensional coupling temis, which yields the following simplification: 


2B 


e, = 


16 


1 A 66 D ll h 


B 


11 


B 


A 11 A 


66 
66 , 


(C30) 
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that is applied to equation (C29). For the subclass of (balanced) antisymmetric cross-ply laminates, 
B12 ~B i6 — B26 — ^66 — 0 , B 2 2 =-B n , and D 16 = D 2 6=0 in addition to the shear-extensional cou- 
pling temis. For this special case, ^ =0 and 


1 - 1/2 


* 


0 



A 12 B ll 


A ll D ll( A ll A 22 _ A 12) e 0 


(C31) 


where 


B 


11 


e n = l- 

A n D n 


(C32) 


For the subclass of balanced, antisymmetric angle-ply laminates, 6^=6 12 = B 22 = 656 = 0 and 
Di6 = D 2 6 = 0 in addition to the shear-extensional coupling terms. For this special case, 


where 


C ° A^D 


(C33) 


66^11 


which agrees with the corresponding equations given by Reuter (ref. 4), and 


<?i = 


B 161 A 1 1 B 26 A 1 2 B 1 6 


2A 66 h \/ A 11 D 11 


A 11 A 22“ A 12 


\ e 0 


(C34) 


Further simplifications can be made to equations (C31) and (C32) for [0/90/... /90] antisymmetric- 
cross-ply-laminate shell walls with an even number of layers that have identical material properties. For 
these laminates, the plies are specially ortho tropic, and their principal material directions are oriented at 
0° and 90° to the cylinder axes in an alternating manner. In particular, the major principal axes of the 
odd-numbered and even-numbered plies are aligned with the x- and 0-axis, respectively, with the first 
ply in the stacking sequence located at the inner surface of the cylinder. Moreover, all odd-numbered 
plies have the same thickness, and all even-numbered plies have the same thickness, but these 
two thicknesses are, in general, different. The laminate stiffnesses are given in reference 18 (see 
pp. 224—226) in temis of the number of layers N, the thickness ratio M, the ratio of the principal elastic 


moduli F = 


% 

Ei 


(for which 0 < F < 1), and the reduced, plane-stress lamina stiffnesses. The thickness 


ratio is defined by 


N-l 

M = X tgc)-*- 

k=l,3,... 


N 


1 


t(k)- 


l (2) 


(C35) 
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where t^) denotes the thickness of the kth ply and 


N-l 


N 


X tfk) + 


2 t (k)-f( t (l) + t (2)) 


(C36) 


k=l,3,... k=2,4,... 


is the total laminate thickness. For the antisymmetric cross-ply laminates, and t^) are the thicknesses 
of the 0° and 90° layers, respectively. Substituting the nonzero laminate stiffness expressions for this 
class of antisymmetric cross-ply laminates that are given in reference 18 into equations (21), (C31), 
and (C32) yields 


0 = < 


[1 -(1-F)Q](M + 1)(M+F) 
(M + F)(l +MF)-[(1 +M)Fv 12 ] 1 2 


1/4 


(C37) 


^0 ~ \/ e 0 


1 + 


\ - 1/2 


/12MF(1-F)v 12 

N^/ [1 -(1 -F)Q](M+ 1)(M+F){(M+F)(1 +MF) [(M + l)Fv 12 ] 2 }e c 


(C38) 


e n = l 

0 ~ t 2, 


12M 2 (1 -F) 2 


N (M + 1) (M + F)[l - (1 - F)Q] 


(C39) 


where 


1 | 8M(M- 1) 

1 +M N 2 (M+1) 3 


(C40) 


and V 12 is the major Poisson’s ratio. For the special, but practical, case of regular antisymmetric cross- 
ply lamination, all plies have the same thickness and equations (C37)-(C39) reduce to 


0 = 


'* v 12 ) 2 


1+F 


-1/4 


(C41) 


^0 _ \/ e 0 ' 1 + 


n/12 v 12 F(1-F) 


- 1/2 


N(1 + F)* / [(1 + F) 2 - (2Fv 12 ) 2 \e ( 


(C42) 




1 -F' 

1 +FI 


(C43) 
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Appendix D 


Example 

Numerical calculations of the nondimensional orthotropy and anisotropy parameters and the bend- 
ing boundary-layer decay length are presented in this appendix for a [(±45)2/304/(0/90)2 ]t unbalanced, 
unsymmetrically laminated cylinder with a radius R = 4 in. For this laminate, the ±45° plies and the 0° 
and 90° plies are made of the Kevlar 49-epoxy material given in table 1 . The other plies are made of the 
AS4/3502 graphite-epoxy material given in table 1. All ply thicknesses are 0.005 in., and the total wall 
thickness is given by h = 0.06 in. 

The laminate stiffnesses were calculated from the formulas given on page 198 of reference 18 and 
are given in matrix fomi by 


A 11 

A 12 

A 16 


4.1678 

1.3051 

1.0787 


A 22 

A 26 

ZZ 


2.4682 

0.3932 

Symmetric 


A 66 


Symmetric 


1.4022 


X 10 5 lb/in. 


Bn 

B 12 

Bl6 


0.7475 

-1.0051 

-0.1288 


B 22 

B 26 

= 


1.2628 

-0.1288 

Symmetric 


B 66 


Symmetric 


-1.0051 


D 11 D 12 D 16 7.8844 2.8842 

D 22 D 26 = 9.3791 

Symmetric D 66 Symmetric 

The numerical value of the nondimensional orthotropy parameter 0 is obtained by substituting the 
appropriate laminates stiffnesses into equation (20a), which gives 0 - 1.0628. To calculate the nondi- 
mensional anisotropy parameter the inverse of the matrix with the subscripted A-temis is needed. 
This inverse matrix is given by 

a ll a i2 a i6 3.4344 -1.4604 

a 22 a 26 = 4.8619 

Symmetric a 66 Symmetric 

The expressions for a 12 , a 2 2 , and a 2 $ are given by equations (C3)-(C5), respectively. Substituting the 
numerical values for B 11? B 12 , B 16 , a 12 , a 2 2 , and a 2 $ into equation (C6) gives b 2 i = 0.0059 in. Similarly, 
substituting the laminate stiffnesses into equation (C7) gives e 0 = 0.9713. Next, the values for the lami- 
nate stiffnesses, a 2 2 , Cq, and b 2 i are substituted into equation (C2) to get s? 0 = 1.1940. 


-2.2325 

-0.2400 xl0“ 6 in/lb 

8.9163 


0.8749 

0.6464 xlO 1 in-lb 
3.0057 
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The numerical value for the first-order correction factor S\ is obtained by first substituting the 
appropriate laminate stiffnesses into equation (C9) to get C\ - 0.0075. Next, the appropriate laminate 
stiffnesses are substituted into equations (C14)-(C22); then, the resulting values are substituted 
into equations (C10)-(C12) to get a 112 = 6.2541 X 10 -10 in/lb, a 122 = -1.5760 X 10 -9 in/lb, and 
a 126 = -5.2861 X 1 0 9 in/lb. Substituting these and the other required values into equation (C13) gives 
bi21 = 3.2255 X 10 -5 in. Similarly, equation (C8) gives C\ - 0.0022. 

The first-order approximation of the anisotropy factor A is obtained by substituting the numerical 
h h_ 

values for p, — , Aq, and into equation (Cl). For the present example, — =0.015. Substituting the 

lv rv 

numerical values for the last three of these quantities and p . = 1.5, which corresponds to the use of the 

Sanders-Koiter equations, into equation (Cl) gives A - Aq = 1.1940. To compute the exact value 

h 

of anisotropy factor A that is given by equation (20b), the laminate stiffnesses, p, and — are 

_ K 

substituted in equations (B31)-(B34) to obtain A 16 = 1.0782 X 10 5 lb/in., A 26 = 0.3927 x 10 5 lb/in., 
A 66 = 1.3947 x 10 5 lb/in., and B 16 = —1.2555 x 10 2 lb. Next, the value for e is computed by 
substituting these results and the appropriate laminate stiffness values into equation (B49); this action 
gives e — 0.9714. Then, the quantities defined by equations (B45)-(B48) are calculated to get a 12 = 
-1.4604 X 10 -6 inAb, a 22 = 4.8619 x 10“ 6 in/lb, a 26 = -0.2401 x 10“ 6 inAb, and b 21 = 0.0059 in. 
Finally, substituting the appropriate numerical values into equation (20b) gives the exact value of 
anisotropy factor A— A o = 1.1940, which is identical with the first-order approximation for A. 


The nondimensional bending boundary-layer decay length, with anisotropy neglected, is obtained 
by specifying a value for the tolerance parameter 8 and then substituting that value and the numerical 
value for the nondimensional orthotropy parameter 0 into equation (19). For 8 = 0.1, this step gives 


d^ 

v/Rh 


- 1.8595. This value corresponds to 


d° 

— =0.2277, which indicates that the bending action attenu- 


ates at a distance of approximately 23 percent of the cylinder radius. The nondimensional bending 
boundary-layer decay length, with anisotropy included, is obtained by substituting the numerical values 


d® ^ 

for — f and A into equation (18). This step gives , = =2.2202, which corresponds to — =0.2719, 
7 Rh sj Rh R 


which indicates that the bending action attenuates at a distance of approximately 27 percent of the cylin- 
der radius. Overall, these results indicate a very small effect of the laminate anisotropy. 
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For all laminates considered, the results show that the differences between results obtained with the Sanders-Koiter 
shell equations, the Love-Kirchhoff shell equations, and Donnell’s equations are negligible. The results also show 
that the effect of anisotropy in the form of coupling between pure bending and twisting has a negligible effect on 
the size of the bending boundary-layer decay length of the balanced, symmetrically laminated cylinders considered. 
Moreover, the results show that coupling between the various types of shell anisotropies has a negligible effect on 
the calculation of the bending boundary-layer decay length in most cases. The results also show that, in some cases, 
neglecting the shell anisotropy results in underestimating the bending boundary -layer decay length and, in other 
cases, results in an overestimation. 
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